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The energy spectrum and positive excess of the hard component of cosmic rays have been 
investigated both in Rome (50 m above sea level) and at Pian Rosa (3500 m above sea level) 
by means of a counter system with deflecting magnetized cores. The results of several measure- 
ments are discussed, and the following interpretation is suggested: (1) A positive excess of 
the order of 20 percent is found in the hard component, in agreement with the results of other 
investigators. (2) The hypothesis of the existence of several types of mesons is not confirmed 
in the lower atmosphere (between 0 and 5000 m). (3) Assuming for the ratio r/yc* between the 
meson proper lifetime + and the rest energy yuc* the value r/yuc* = 3 X 10-* sec./Mev, one has to 
assume that the energy spectrum of mesons at the top of the atmosphere does not follow a 
power law E~*-* for low values of the energy (E <4X 10° ev). 


1. INTRODUCTION 


N the years 1941-43 we have performed in the 

Physical Institute of the University of Rome 

a number of experiments on the mesonic com- 

ponent of cosmic rays by means of magnetized 

iron cores. We present here an outline of this 

work: some details that will be omitted may be 
found elsewhere.! 

The first attempts to detect the deflection of 
cosmic particles by the magnetic field in an iron 
core were made in 1930 by Rossi? and Mott- 
Smith* independently. Their apparatus were 


* A fund for researches on cosmic rays was granted by 
the Consiglio Nazionale delle Ricerche. 

1G. Bernardini and M. Conversi, La Ric. Scient. 11, 858 
(1940). G. Bernardini, M. Conversi, E. Pancini, and G. C. 
Wick, La Ric. Scient. 12, 1227 (1941). M. Conversi and 
E. Scrocco, N. Cim. I, 372 (1943); Bernardini, Conversi, 
Pancini, and Wick, Phys. Rev. 60, 535 (1941). 

* B. Rossi, Nature 128, 300 (1931). 

+L. M. Mott-Smith, Phys. Rev. 35, 1125 (1930); 37, 
1001 (1931); 39, 403 (1932). 


rather different, but similar in two respects: 
both used a counter telescope to detect the par- 
ticles, and in both cases the effect investigated 
was essentially a differential one, depending on 
the charged particles in the radiation being of 
one sign only or, at least, prevalently so. 
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Fic. 1. Arrangement of counters and magnetic cores. 
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Fic. 2. Coils of pair of cores (A’A” or B’B’’). 


Mott-Smith obtained a negative result.‘ Rossi 
on the other hand, by means of a telescope con- 
sisting of only two counters found an indication 
of a small effect, corresponding to an excess of 
the positive over the negative particles; the 
evidence, however, was partly contradictory and 
was not regarded as definite by the author. 
Looking at these results in the light of later 
knowledge we can see that: (a) the hypothesis 
that the effective field inside iron is not given 
by the induction vector B and is considerably 
smaller than B (this, of course, would easily 
explain the negative result of the deflection ex- 
periments of Mott-Smith and Rossi) has been dis- 
missed by a theoretical discussion by v. Weiz- 
sicker; (b) an alternative explanation is offered 
by Blackett’s result,* namely that the radiation 
consists of positive and negative particles in 
approximately equal numbers; (c) while this is 
probably the correct explanation, we must not 
forget that there is a certain discrepancy between 
Blackett’s statement and the subsequent investi- 
gations of Leprince-Ringuet and Crussard,’ 
Jones,* and Hughes.® Especially the last two 
authors have given quite definite evidence of the 
existence of a small but real excess in the num- 
ber of the positive over the negative mesons; 
the excess is of the order of 20 percent, and 
according to this, one should find in an experi- 
ment like that of Rossi a small but definite effect. 
It would seem, therefore, that the subject 
deserves further investigation; but there are 
more cogent reasons too. The positive excess in 
the meson spectrum is probably connected with 
the positive nature of the primary radiation 

4In a su uent investigation, however, Danforth and 
Swann, Phys. Rev. 49, 582 (1936) have found a small but 
real effect by means of an apparatus not very different 
from that of Mott-Smith. 

5 C. F. Weizsacker, Ann. d. Physik 17, 869 (1933). 

*P. S. M. Blackett, Proc. Roy. Soc. 159, 1 (1937). 
(1937) Leprince Ringuet and J. Crussard, J. de. Phys. 8, 207 


5H. {tie Rev. Mod. Phys. 11, 235 (1939). 
®D. J. Hughes, Phys. Rev. 57, 592 (1940). 


that is revealed by the azimuthal effect; a 
study of the variation of the positive excess with 
height would probably be interesting and might 
throw some light on the process of creation of 
the mesons. Now the Wilson chamber method 


a) b) 


Fic. 3. Qualitative shape of trajectories: (a) “‘cc’’ arrange- 
ment; (b) “cd” arrangement. 


1 T. H. Johnson, Rev. Mod. Phys. 11, 208 (1939). 
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employed by all the authors mentioned above 
is obviously unsuitable for this purpose owing 
to the rather heavy equipment it requires. 

Bearing this in mind, we have tried to improve 
Rossi’s experiment, and we have planned an 
apparatus consisting of a triple coincidence 
system with four iron cores; we have found that 
it gives a quite measurable effect, depending 
upon the positive excess. We hope to use it for 
an investigation of the positive excess in the 
upper atmosphere. 

This apparatus also allows an entirely different 
use since the magnetization of the cores can be 
arranged in such a way as to cause a deflection 
effect which is not differential. That is, the effect 
depends on the energy of the particles, but not 
on their sign; accordingly, we have found in this 
case a substantially larger effect, which can be 
used for investigation of the energy spectrum of 
the particles." 


2. APPARATUS 


Our apparatus is shown schematically in Fig. 
1; it consists primarily of three coincident coun- 
ters having a radius r=1 cm and useful length of 
about 50 cm. The distance d between the axis 
of the intermediate counter ‘‘2’’ and that of 
either counter “1” or “3’’ is 22 cm. The two 
pairs of iron cores A’, A” and B’, B” are magne- 
tized by coils as is shown by Fig. 2. The field 
in the cores is parallel to the axis of the counters 
and has opposite directions in A’, A’’ (similarly 
in B’ B”’). The magnetic circuit is closed by 
iron bars applied at both ends of A’, A” (or 
B’, B’’). Thus each pair of iron cores, for instance 
A’, A”, acts, roughly speaking, like a cylindrical 
magnetic lens which concentrates the particles 
of one sign. We shall hence-forward speak of 
this ‘“‘lens’’ as convergent (or briefly ‘‘c’’) or 
divergent (briefly “‘d’’) for positive (c+ or d+) 
or negative (c— or d—) particles. 

"The main results of the present paper have already 
been indicated in a preliminary communication (Phys. 
Rev. reference 1). After the present paper was completed 
we have learned from a private communication that 
F. Rasetti has recently succeeded in obtaining a deflection 
effect by an arrangement of the type used by Mott-Smith. 
His results are in accord with the theoretical prediction 
that the effective field inside iron is B and with the ex- 
istence of a positive excess. All this seems to be in the best 


ge ei with our own previous conclusions and with 
of the present paper. 


Fic. 4. Modification impressed u the total en 
spectrum (positive+negative mesons) by the action of the 
instrument with either arrangement “‘cc”’ or “‘cd."” 


Now the upper (A’, A’’) and lower (B’, B’’) 
pairs can be magnetized in the same way or not, 
so that on the whole we have the following four 
possibilities: 

cc+ (=dd—); dd+(=cc—); 
cd+ (=dc—); dc+ (=cd—). 


The two alternatives cc and cd are sche- 
matically shown in Fig. 3. 

With the “cc” arrangement, the number of 
positive particles detected by the counter tele- 
scope is increased, the number of the negative 
particles is decreased, and the reverse takes place 
with the ‘dd’ arrangement. The difference 
Nect—Nee— between the number of triple co- 
incidences, N..+ or N..— registered in a given 
time with the two arrangements depends on the 
asymmetry of the meson spectrum; i# may be 
used to measure the positive excess, as a more 
quantitative analysis will show. 

We shall later develop a theory of the action 
of the instrument. Figure 4 gives an idea of the 
modification impressed upon the total energy 
spectrum (positives+negatives) by the action 
of the magnetic lenses with either arrangement 
cc or cd. In this case the positive excess has been 
neglected; the broken line in Fig. 4 shows the 
energy spectrum at sea level in the vertical 
direction and in open air. This spectrum has 
been evaluated by assuming, as usual, a power 
spectrum for the mesons at the creation level 
in the upper atmosphere, and by taking into 
account the disintegration of the mesons, and 
by assuming for the mean lifetime of +r of a 
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TABLE I. Measurements of the positive excess. 


SCROCCO, AND WICK 


Meters 
above sea level 


Threefold coinc./hour 


Nee+ 


Double coinc./hour 


c+ 


Ne- 


50 
(Roma) 


81.141.2 
83.1+1.2 
88.5+0.9 
97.5+1.0 
103.3+1.5 


170.0+2.1 
180.5+2.5 
190.5+2.1 
221.0+41.8 
223.042.7 


77.741.2 
77.2+1.1 
80.5+0.9 
84.7+0.9 
91.2+1.3 


175.0+2.2 
172.0+2.5 
174.5+2.0 
191.0+2.0 
200.0+2.9 


235.4+42.4 
243.6+2.4 
245.2+1.9 
253.7+2.0 
260.2+2.8 


615+8.0 
597+9.5 
609 +4.5 
664+4.0 
678+7.5 


234.442.3 
235.7+2.3 
236.2+1.8 
242.341.9 
250.5+2.5 


601+8.0 
589+9.5 
583+4.5 
624+4.0 
63147.5 


— 


. 


NOP 


meson a value:” r/yc?=2X10-® sec./Mev (yu is 
the meson mass) with uc? = 100 Mev. 

The full lines give the spectrum as modified 
by the field (and the absorption in the iron 
cores) when the field B in the two pairs of iron 
coils is B= 17,000 gauss. It can be seen that the 
spectrum is noticeably modified in the energy 
interval, between ev—the minimum 
energy for a meson to get through the iron bars— 
and 410° ev. The spectrum becomes richer in 
low energy particles with the cc arrangement; 
the reverse takes place with the cd, or dc, arrange- 


ment. 
One can see, also that the energy range of 
field sensitive particles is about the same as for 


the Wilson chamber technique." 


12 This value of r/yc* is, in our opinion, rather too low; 
-but we have chosen it in our calculation because it is 
intermediate between the several values found by different 
experiments. 

43In comparing the sensitivity of our apparatus with 
that attained in the Wilson chamber, it is well to remember 
that in the latter case the limit to the measurement of 
large energies is not set by the difficulty of measuring the 
small curvatures involved, but rather by the spurious 
curvatures produced by irregular motions of the chamber 
gas. In our case this cause of error does not exist; unfor- 
ge it is superseded to a certain extent by a similar 
effect due to scattering of the mesons in the iron cores. 
This causes a spurious curvature of the trajectories, and 
thus tends to obscure the effect of the magnetic field. 
The following evaluation shows, however, that the change 
in the action of the magnetic “‘lens’’ by the scattering, 
must not be very large (and, at any rate, the spurious 
curvatures due to scattering vary as the reciprocal of the 
particle energy; that is, they become smaller at the same 
rate as the true curvatures, which is not the case in the 
Wilson chamber !). 

The order of magnitude of the spurious curvatures 
introduced by scattering can be evaluated by means of 
Williams’ formulae, Proc. Roy. Soc. 169, 531 (1939), which 
give the mean deflection 3 suffered by a particle in going 
through a plate of thickness ¢, of atomic number Z, 
density p, atomic weight A. Consider, for instance, a 
particle moving in a plane normal to the axis of the counters 
and call @ the component of the deflection parallel to this 
plane—this is in fact the component in which we are 


The Wilson chamber, of course, allows a much 
more detailed analysis of the energy spectrum 
to be made, but, this advantage is attained only 
by the much greater labor which the method 
requires. 

On the whole we think that the two methods 
are not mutually exclusive, but that they can 
usefully integrate each other. , 


3. RESULTS 
(a) The Positive Excess 


These measurements were made partly in the 
Physical Institute of the University of Rome 
and partly at Pian Rosa (Cervinia) 3500 m 
above sea level in a wooden hut." 

During the measurements the two arrange- 
ments cc+ and cc— were, of course, systemat- 
ically alternated; also, several values of the 
magnetic field were tested. The electric pulses of 
counters were sent to a coincidence system which 
registered besides the triple coincidences also the 
double coincidences between the two upper 
counters. 
mainly interested, as it is indistinguishable from a deflec- 
aed to the magnetic field—we find, according to 


600Ze(Nt)! 
p= [3.7+0.28 log | 


’ where E is the energy of the particle in e-volts, 8c its 


velocity (8~1), and N is the number of atoms per cm’, 
= iron we find: 3=~((1/t)/E)10’, when ¢~20 up to 


cm. 

The deflection caused by the magnetic field is: #3’ 
=300Bt/E, which becomes, when B=15000 gauss: 
0’=(t/E)10°. Putting cm, one can see that is 
about three times as large as #, independently of energy. 

4 A series of measurements in a plane flying at 9000 m 
above sea level was also begun, but unfortunately a flying 
accident destroyed the apparatus before a sufficient set 
of data had been secured. 
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These could be compared with Rossi’s results, 
although the distance between the counters and 
the counter diameter are not quite the same. 
Our results are shown in Table I. 

The effect 6 in the sixth column is defined by: 


Nee+— Nec— 
Nees + 


Il 
N 


(1) 


The effect 6’ in the seventh column is defined in 
a similar way for the double coincidences. 

The value of 6 has been evaluated without 
taking into account the effect of showers, which 
had been shown in the usual way (that is by 
shifting the upper counter horizontally, out of 
the solid angle defined by the remaining two 
counters) to be negligible. 

The correction for random coincidences is 


"negligible in the case of the triple ones; the double 


coincidences due to chance were about 40/hr. 
at Pian Rosa as well as in Rome" and they are 
already taken account of in the N+ and N— 
values given in Table I. 

In Figs. 5 and 6, the values of 6 and 4’, found 
in Rome and at Pian Rosa, are plotted against 
the field B. The effect on double coincidences, 8’, 
is seen to be quite a bit smaller than on 5; this is 


| | 


Fic. 5. Results of the measurements of the positive excess 
performed in Rome (50 m above sea level). 


“The building material of the Institute is weakly Fic. 6. Results of the measurements of the positive excess 


radioactive and increases the natural pulses of counters. 
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understandable since doubling the distance be- 
tween the extreme counters not only increases the 
angular definition of the beam selected by the 
telescope, but also doubles the deflection caused 
by the field. In fact, the theory of the apparatus 
(Section 4) shows that the sensitivity increases, 
roughly, with the square of the distance between 
the counters. 

Whatever the reason may be, the effect 8’ 
which we have found with the double coinci- 
dences is much less reliable than 6. 

Considering the values of 6 found at Rome and 
Pian Rosa, one can see that the results are quite 
similar, indicating that the positive excess in the 
meson spectrum does not materially change 
between sea level and 3500 meters. 

With regard to the dependence of 6 on B, 
one may remark that 6 can never be larger than 
the positive excess e¢, if ¢ is distributed uniformly 
over the energy spectrum. Indeed one can see 
that 6 would finally become equal to e« if the 
action of the “lenses’’ would become strong 
enough to exclude the particles of one sign 
completely. The 6(B) curve should therefore 
resemble a saturation curve. But with the field 
values we have reached, there is no indication 


of a bend. Our results, however, are not incon- 


sistent with the hypothesis that for higher values 
of B, 6 might tend to a saturation value. 
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performed at Pian Rosa (3500 m above sea level). 
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TABLE II. Measurements on the meson energy spectrum (B = 17,000 gauss). 
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Absorbing 


Threefold coincidences 


Meters layer hour 
Series above sea level (g/cm?) Arrangement (Neca; Neda) at a 4=a4—a_ 
cc+ 104.0+1.4 
50.8+3.0 
cd+ 69.6+1.0 
1st 50 0 42.342.0 16.1+3.9 
(Roma) cc— 94.6+1.2 


2nd 50 250 cd+ 


70.7+0.9 


91.941.4 
62.7+1.0 
80.0+1.4 
61.0+1.0 


233.743.8 


48.0+3.3 


40.0+2.3 15.8+4.6 


32.242.2 


67.244.5 


cd+ 145.6+2.7 


(b) Measurements of the Meson 
Energy Spectrum 


The results of these measurements, Table I], 
are grouped in three separate series, each series 
comprising measurements made at the same 
height above sea level and with the same ab- 
sorber thickness above the apparatus. 

In the first two series of measurements, made 
in Rome, the direction of the current was 
systematically changed in the coils around A’A”’ 
~ and independently in the coils around B’B’’, so 
that all four arrangements cc+, cc—, cd+, cd— 
were alternated. 

With the mean values 


(2) 
the effect we were looking for was defined by 
a =100(Nec— Nea) / Nea: (3) 


This is given in the next to the last column of 
Table II. In the third series (at Pian Rosa 3500 
m above sea level) the current in the A’A” coils 
was never changed so that the upper ‘‘lens’’ was 
always convergent for positive particles. The 
effect was, therefore, defined by 


ay =100( — Neay)/Nea+ (4) 


The second series of measurements was planned 


after Nielsen, Ryerson, Nordheim, and Morgan" 
had announced that a new careful investigation 
of the proper lifetime of a meson had given a 
value 1.25 X sec./Mev, considerably 
lower than all values previously found. 

With this value one can see that the maximum 
in the energy spectrum of the mesons arriving 
at sea level in the vertical direction should lie 


above the minimum energy of the mesons which 


are not stopped in the 40 cm Fe of our apparatus. 
Now if a thick absorber is placed above the 
counter telescope, it will cut off the lower part 
of the energy spectrum; for a properly chosen 
thickness of the absorber one can see that the 
resulting spectrum—that is the spectrum im- 
pinging on the apparatus as modified by the 
absorber—is appreciably softer than the un- 
modified spectrum. 

We might then expect the effect a to be in- 
creased; this is also confirmed by the more 
detailed considerations of Section 5. Accordingly 
we built a brick wall, under which our apparatus 
could be placed, covering the solid angle defined 
by the counter telescope. 

The vertical thickness of the wall was 250 g/ 
cm’, and the results are grouped in the second 
series. 


The effect of showers has been deducted from 


16 W. M. Nielsen, C. M. Ryerson, L. W. Nordheim, and 
K. Z. Morgan, Phys. Rev. 59, 547 (1941). 
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the results included in the first two columns of 


Table II; the correction was, however, very 
small. The data given in the last column can be 
used for a further evaluation of the positive 
excess. 


4. THEORY OF THE INSTRUMENT 


We shall now try to draw from the experi- 
mental data a quantitative evaluation of the 
positive excess as well as some more detailed 
information about the meson spectrum. 

The direction of motion of an incoming 
particle shall be specified by the polar angles 3 
(colatitude) and ¢ (longitude). As polar axis 
we chose the axis of a counter, or the direction 
of the magnetic field inside the cores. Let further 
y=0 indicate the vertical plane containing the 
three counter axes. Put 


dn=J(E, 8, ¢) sin ddddgdE, 


for the number of incoming particles per second 
and per cm?, of energy E within the interval dE, 
and of direction 3, g within the solid angle 
sin dddd¢. 

For particles of a given sign, there must be a 
“cross section” of the instrument =(E, 8, ¢) 
such that the number dN of incoming particles 
within the specified intervals of energy and direc- 
tion, that cross the useful volume of the three 
counters, is given by 


dN=dn (E, ¢). 


Naturally = depends also on the field B, on the 
geometrical constants of the instrument, and on 
the field “‘arrangement” (whether cc, or cd, etc.). 
Since the instrument is sensitive (or +0) 
only to particles with | ¢|«1 we may neglect 
the dependence of J on ¢ and use for the total 
number of mesons crossing the counter system 


N= f f J(E, 8)dE sin odo f (E,9, (5) 


Now a more detailed investigation’? shows that 
the second integral can be put in the form 


f (6 


™ Bernardini, Conversi, Pancini, and Wick, reference 1. 


where o and y are dimensionless measures for 
the energy E and the stopping power L of iron, 
respectively, 


E; v=—L. (7) 


Here d is the distance between the axis of the 
central counter and either the lower or the upper 
counter, 7 is the radius of a counter, e is the 
charge of the particle. In Eq. (6) y: is a simple 
function which needs not to be known exactly; 
the function f(c, 7) shall be indicated as the 
“characteristic function’’ of the instrument and 
expresses the dependence of the cross section on 
all parameters involved through the relation (7). 

We shall now suppose that the dependence of 
the spectrum J(E, #) on energy and # can be 
separated so that we can put 


J(E, 8) (8) 


This hypothesis can be justified approximately; 
remember that the intensity of mesons of energy 
E, incident under a zenith angle ¢, is well repre- 
sented! by the formula 


m n 
[z +7] dE, (9) 
cos 


with a proper choice of the constants m and n. 
Now since in our measurements the counter 
telescope was vertical, and remembering that 
| <1, we have: cos {sin 3 and putting 


d*Be Ld 
a= =—; y.(8)=sin"d; 
r Y 
=(ao+m)-", (10) 


we really find formula (8). 
Substituting (6) and (8) into (5) and integrat- 
ing over #, we find 


N= f &(a)f(c, (11) 


neglecting a constant factor which is unessential 
in the following development. 

We have termed f the characteristic function, 
but in fact, there are four such functions, corre- 


18 Conversi and Scrocco, reference 1. 
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TABLE III. Results of the a and a, calculated effects. 


tat” beh 800” bat 
1 38.7 17.9 44.8 24.4 56.0 36.2 


2 79.5 41.8 78.8 38.6 119.0 66.8 
3 100.0 51.1 93.8 47.5 152.1 81.5 
4 111.7 55.6 102.6 51.3 171.1 92.5 


sponding to the four possible arrangements, cc, 
dd, cd, dc (“‘convergent” ‘‘divergent’’ referring 
of course to particles of the specified sign). We 


may indicate these functions by fee, faa, *** 


respectively; if they are known, we may proceed 
to a theoretical evaluation of (1), (3), and (4). 
How these functions can be determined has 
been discussed elsewhere.'7:!§ We shall only recall 
some results. 
We shall put 


(12) 


which is the value of the ‘‘energy”’ « possessed 
by a particle when it reaches the central counter, 
that is after suffering one-half of the total energy 
loss in crossing the iron cores. The theory then 
shows that one has, with sufficient approxima- 
tion 
See(o, ¥) 0); faa(o, =faa(é, 9) ; 
3Lfea(o, +faclo, J=fealé, 0) =fac(é, 0). (13) 


These formulae have a rather simple meaning; 
the first two mean, for instance, that the char- 
acteristic function can be evaluated by neglect- 
ing the energy loss (that is putting y=0 in 
f(c, y)) provided the energy ¢@ instead of the 
initial energy o is taken. Then f(c, y) may be in- 
dicated more briefly with f(¢). We have found 


feel) = 1+ (1/6) — ; 
1—(1/6)+(1/46*); (#>3) 


[¢<}; B> ~4.10° gauss] 


1+¢ 1 G 
fea(4) =fac(é) = + 262 Em 4g? 


where &,, is amongst the three real roots of 
the one which satisfies —1<&, <0. 
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The value of the coefficient k must still be 
considered. The term containing k has been 
introduced to take into account in a rough way 
the effect of scattering within the iron cores, 
which could not be included in the theory. It is 
natural to assume that this effect will be im- 
portant especially at the low energies and espe- 
cially in the ‘‘cc’’ case when the cross section of 
the counter system tends to increase to very 
large values with decreasing energy. In our 
simplified theory this increase continues until 
the cross section is suddenly brought down to 
zero, because the particles are stopped in the iron 
before reaching the lowest counter. Scattering, 
however, must set a limit to this increase before 
the stopping limit is reached, and this effect 
may be represented by a term of the type con- 
sidered. The choice of the value of & has no 
considerable influence on the evaluation of the 
positive excess (see Section 5(a)); unfortunately 
this is not so in the case of the meson spectrum. 
We have avoided the difficulty by the assumption 
that the proper value of k can be determined by 
comparing our experimental results in air and 
at sea level with calculations based on formulae 
(11) and assuming as energy spectrum one that 
can be derived from the absorption curve in lead 
measured by Clay and van Gemert. For further 
details see one of our previous papers;'® the 
result is: k=0.13+0.01. 


Fic. 7. Theoretical 5(B) curves (see formula (1)) under 
the assumptions of a uniformly distributed positive excess 
(curves k=0 and k=7) and of an excess localized about 
E’=3X 10° ev. 


th 
A 
U 
Ca 
we 
wl 
(F 
$= 
wi 
0.7 
ass 
ki 
48 
vei 
an 
me 
AX 
| 
05 hy 
the 
to 
bil 
thi 
Nes 
are 


MAGNETIC DEFLECTION OF COSMIC RAYS «117 


5. DISCUSSION 
(a) The Positive Excess 


Remembering the definition of N.., Naa, by 
means of Eq. (11) we have 


Noa = Ido, 
(14) 


the meaning of the symbols being self-evident. 
Assume the positive excess ¢ to be distributed 
uniformly over the spectrum 


( 
(15) 
= i( 1 a). 


Using the third of Eqs. (10), with »=%3, we have 
calculated N.., Naa as functions of B. From these 
we derive the 6(B) curve represented in Fig. 7,'° 
which is in good agreement with the approxi- 
mately linear increase of the experimental 6 
(Figs. 5, 6). We have found experimentally 
§=0.15 for B=17,000 gauss; comparing this 
with the two curves of Fig. 7, we find e=0.19 or 
0.23 according as we make for k the two extreme 
assumptions k=0, and k=. The exact value of 
k is therefore rather unimportant, and we may 
conclude that e=20 percent approximately in 
very good agreement with Hughes and Jones. 

The agreement between the result of Hughes 
and Jones and our own gives a plausible argu- 
ment in favor of B as the effective field inside 
iron. 

Subsequently we have also examined different 
hypotheses about the distribution of the positive 
excess over the spectrum, but we have found that 
the shape of the 6 vs. B curve is rather insensitive 
to such changes, so that closer investigation of 
the shape of this curve does not open the possi- 
bility of investigating the distribution of the 
positive excess over the spectrum. The only 
thing we can say is, that if the excess were 


? We have not included in the figure the N,. vs. B and 
* Naa vs. B curves, as derived theoretically; but they also 
are in good accord with the data of Table I. 


TABLE IV. Results of the comparison of theoretical curves 
for a(r/yc*) and experimental results. 


Series 
+0.40 
+0. 

2nd 2.25 0.38 
+0. 

3rd 2.10; 5°35 


located in the low energy part of the spectrum, 
say around E<~10° ev the value of 6 to be 
expected (for a positive excess of 20 percent) 
would be considerably larger than observed. 


(b) Energy Spectrum 


In this case, we assume, for greater accuracy, 
the meson spectrum to be given by the product 
of (9) by the probability W (E, h) for a meson of 
energy E to reach the depth A without dis- 
integrating. In this way the effects a and a,, 
defined by (3) and (4) are seen to depend on the 
7/yc? ratio, which appears in the expression for 
W (E,h). If j is the stopping power of air, we 
shall have m=hj. Further we have chosen 
n=2.87 as the most accurate value known to 
us.2° For W (£,h) the usual expression was 
used, and by means of (7) the spectrum was 
expressed as a function of the dimensionless 
energy o. 

Assume the positive excess to be uniformly 
distributed over the spectrum. Writing the ex- 
plicit expressions of N.g and Na, (analogous to 
(20)) and putting for shortness sake: 


F(a) = fee(o) +faalo) J; 
we shall have (see Eq. (2)), 


f Fle)de: 
(16) 


Ne= f 


while (4) can be expressed approximately” in 
20 A. Emhert, Zeits. f. Physik 106, 751 (1937). 


That is neglecting a term [fae(o, y)—fealo, y) ldo 
compared to Nea. 
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Fic. 8. Comparison of theoretical a(r/uc*) curves (see 
formulas (2)—(4)) and experimental results. 


the form 


17 
f (o)g(e)de. (17) 


ed 


The integrals in (16) and (17) have been 
evaluated numerically, for each of the three series 
of measurements, for two values of k (k=0 and 
k=1/8) and for four values of the ratio r/yc? 
(=1; 2; 3; 4-10-® sec./Mev). The results are 
given in Table III. Comparing the calculated 
a(r/uc?) curves (these can be drawn from the 
results of Table III also for the k values con- 
sidered in Section 4, remembering the linear 
dependence of a on k) with the experimental 
results of Table II, we obtain for r/c? the values 
contained in Table IV. The essential hypothesis 
made is that the meson spectrum at the top of 
the atmosphere can be represented by the 
law E-?-87, 

The calculated curves of Fig. 8 confirm the 
statements already made in Section 3(b) about 
the influence on @ of an absorber 250 g/cm? 
thick, if one assumes for 7/yuc? a value as low as 


PANCINI, 


SCROCCO, AND WICK 


that found by Nielsen e¢ a/;"* the effect observed 
experimentally has not the sign to be expected 
on the assumption of a value as low as that. 

The values for r/uc? derived from the three 
series of measurements are in surprising agree- 
ment with each other, especially considering the 
fact that the agreement is not sensitive to a 
change in the initial spectrum. This seems to be 
one further reason to exclude a complex nature 
of the meson component,” at least between sea 
level and 4000 m. 

Let us now briefly discuss our results. The 
evaluation of r/uc* was based on formula (11); 
hence it is subject to errors due to (a) the ap- 
proximations adopted in the theoretical treat- 
ment of the characteristic functions and (b) 
deviations of the actual energy distribution at 
the top of the atmosphere from the distribution 
assumed, formula (9). 

As regards the first cause of uncertainty a 
discussion of the hypotheses involved has been 
given elsewhere'® but it is rather hard to estimate 
accurately the possible error. We feel, however, 
that to assume that this is the source of the 
difference between our present results and the 
value, +/uc?~3X10-° sec./Mev, which is indi- 
cated by the most recent and reliable measure- 
ments by the anomalous absorption method,” 
would be to stretch the limits of error of the 
calculations more than is probable. 

We shall, therefore, examine the second pos- 
sible cause of discrepancy. Obviously one can 
interpret the discrepancy by the assumption, 
that the meson energy spectrum on the low 
energy side is less rich in particles than corre- 
sponds to (9). According to this formula the 
spectrum is given by the product of (9), with 
n= 2.87, by the probability of survival W (E, h). 
This probability was evaluated on the usual 
assumption that the creation of mesons takes 
place at a depth of 100 g/cm? below the top. It 
is well to remark that the value of W (E,h) 
at sea level is not very sensitive to changes in 
the depth of the creation layer. Turning now to 
formula (9), we notice that it is supported by 
the work of Blackett, Jones, Hughes, and others, 
for n=2.9 and for E=values comprised between 


ae Weisz, Phys. Rev. 59, 554 (1941 ). 
% G. Bernardini, Zeits. f. Physik 120, 413 (1943). 
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some 10° ev and the maximum energies measur- 
able by the Wilson chamber method. The inten- 
sity or depth curves of Ehmert and others 
indicate that the validity of the formula extends 
to much higher energies. Since moreover the 
exact shape of the spectrum above 10! ev is 
immaterial for the present purpose, we may 
assume the formula E-*-*’ to represent the initial 
spectrum correctly for all E-values larger than 
a certain E*. Considering that in crossing the 
atmosphere a meson loses about 2.5X10° ev, 
we shall put E*~4X 10° ev. 

For E<E* the E~-*’ law is not sufficiently 
supported by any experimental data. 

There is even some indication to the contrary, 
for the energy spectra of Hughes and Jones seem 
to be less rich on the low energy side than one 
would expect assuming a creation spectrum 
according to (9), this conclusion being practically 
unaffected by a variation of the assumed height 
of the creation layer between 50 and 200 g/cm’. 

Such a conclusion would be in agreement with 
our observations and would not lead to any 
difficulty in connection with present knowledge 
about the meson radiation in the upper atmos- 
phere. As a matter of fact, the experiments of 
Schein, Wollan, and Groetzinger,™* and of Herzog 
and Bostick,”®> have revealed that a large number 
of rather slow mesons are created in the upper 
atmosphere; but these particles are in general 
much too slow to reach our apparatus and be 
recorded by it. 

One further remark is that the meson spectrum 
must probably be affected in some way by the 
cutting off in the spectrum of the primary pro- 
tons (?) by the magnetic fields of the earth and 
sun. This would cause a falling off on the low 
energy side. On the other hand one may com- 
pare these suggestions with the results of the 
measurements of 7+/yc? by the anomalous absorp- 
tion method based on integral effects. These 
require hypotheses on the creation spectrum of 
mesons which are quite analogous to those which 
are necessary to discuss our experimental results. 

As a matter of fact, in the anomalous absorp- 
tion method, one finds that the agreement 
between the r/yuc? values based on differential 

*™M. Shein, E. O. Wollan, and G. Groetzinger, Phys. 
Rev. 58, 1027 (1940). 


*G. Herzog, and W. H. Bostick, Phys. Rev. 58, 278 
(1940); 59, 122 (1941). 


effects, which are independent of any hypothesis 
concerning the spectrum, and those based 
on integral effects, which are not, is in general 
rather poor. This points to some hypothesis 
such as we have indicated. Two notable ex- 
ceptions are given by the measurements dis- 
cussed by Bernardini and Festa** and by the 
ionization chamber measurements of Neher and 
Stever.?’ As regards the first exception, however, 
one may remark that the agreement found may 
be owing to the evaluation method employed 
by those authors, a method whose precise aim 
was to eliminate the uncertainty of the low- 
energy side of the meson spectrum from the 
final result. 

As regards the experiment of Neher and 
Stever, one has to remark that it refers to the 
meson radiation incident from the whole hemi- 
sphere; consequently the mean path in air of 
the particles was much larger than with a ver- 
tical counter set. This naturally increases the 
importance of the high energy side of the spec- 
trum, where no large discrepancies from the 
E-*-87 law need be assumed. 

Let us finally turn our attention to the data 
of the last column in Table II. These allow (see 
also Section 3(b)) a further evaluation of the 
positive excess ¢, on the assumption that it is 
uniformly distributed over the spectrum. In fact, 
putting 


1 
J &(o)g(o)de, 


we get from (17) and from the analogous expres- 
sion for a_ 
A=ue. 


The evaluation gives: 


_ ) (18.344.4) percent in open air; 
’ (18.6+5.4) percent under a layer 250 g/cm? 
thick, 


in good agreement with the e-value derived from 
Table I. There is as yet no evidence that the 
250 g/cm? layer causes an appreciable change in 
the positive excess. 
We finally wish to emphasize the possibilities 
%* G. Bernardini and C. Festa, Atti d. R. Acc. D’Italia 
4 (1943). 


H. V. Neher and H. G. Stever, Phys. Rev. 58, 766 
(1940). 
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which the present method may have in the in- 
vestigation of the meson spectrum at great 
heights (notice in this connection how the slope 
of the a vs. r/uc? curves increases with increasing 
height). An experiment of this kind by Bernardini 
and Scrocco was in progress in 1941, but was inter- 


rupted by a serious flying accident to the airplane 
which carried the apparatus. As a consequence 
the apparatus was gravely damaged, and for 
this and other reasons connected with Italian 
events, the experiment was given up for the 
time being. 
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Disorder Scattering of X-Rays by Local Distortions 


H. EXKsTEIN 
Armour Research Foundation, Chicago, Illinois 


(Received July 16, 1945) 


The disorder scattering (background) owing to local distortions is caused mainly by the 
elastic strain field surrounding the distorted zone rather than by the misfit atoms in the zone 
of distortion. As an example, the scattering caused by a distortion of spherical symmetry is 
calculated. The result is extended for more general types of distortions. The background 
intensity increases strongly in the neighborhood of a line, and it tends toward a constant value 
for small scattering angles. This explains previous observations on rolled copper. 


1 


explain the plastic properties of solid 
bodies, it is necessary to assume local lattice 
distortions or irregularities. Most of the evidence 
for the assumed models of dislocation is indirect. 
Several papers have attempted to connect the 
observed line intensities and widths with the 
internal distortions.’~* Little attention has been 
given to the continuous background scattering 
arising from the distortions. The experimental 
conditions for this type of research are not easy, 
because other effects such as fluorescent radia- 
tion, Compton radiation, diffuse scattering by 
air, coherent temperature scattering and back- 
ground owing to the continuous x-ray spectrum, 
all contribute to the observed background. Yet 
experiments by Guinier’ have shown the effect of 
plastic deformation of copper upon the back- 
ground to be appreciable. If the scattering caused 


a9 acre and H. Mark, Zeits. f. Physik 61, 435 


as $3) W. Brindley and F. W. Spiers, Phil. Mag. 20, 882 
3 U. Dehlinger, Zeits. f. Krist. 65, 615 (1927). 
*W. Boas, Zeits. f. Krist. 97, 354 (1937). 
5 U. Dehlinger and A. Kochendérfer, Zeits. f. Metallkunde 
31, 231 (1939). 
* A. Kochendérfer, Zeits. f. Krist. 101, 149 (1939). 
7 A. Guinier, Comptes rendus 208, 894 (1939). 


by local distortions can be isolated, it will provide 
more useful information than the study of the 
lines; first, because this scattering is entirely 
absent in an ideal crystal, while the line intensity 
decreases only by a few percent when distortion 
is introduced, and secondly, because slowly vari- 
able internal stresses affect the line, but not the 
background. 

The following theoretical discussion gives a 
correlation between local distortions and disorder 
scattering, which explains some of Guinier’s 
results, and may be helpful for further experi- 
mental work. 

For a simple Bravais lattice, the scattering 
intensity is given by: 


I=I,| F|?| exp (2iR-r,) |’, (1) 


where J, is the scattering intensity from one 
electron, F the atomic form factor, 


k and kp the wave vectors of incident and 
scattered wave, and r, the radius vector of the 
nth atom. For r,, we write 


Tn =AntUn, (3) 
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where a, is the radius vector of the mth lattice 
int. 

In physical lattice distortions of the type 
usually assumed, the displacements of the atoms 
from the regular lattice points are of the order of 
magnitude of the atomic distances in the immedi- 
ate neighborhood of the distortion which extends 
over only a few atomic distances. The number of 
atoms within this zone is small, and they cannot 
be expected to give a large contribution to the 
disorder scattering. Outside the distortion zone 
itself, the displacement of atoms is of the elastic 
type, i.e., the displacements are small as com- 
pared to the atomic distance. While the displace- 
ment decreases rapidly from the origin of the 
perturbation, the number of atoms affected is 
large, and we may expect an important x-ray 
effect from them. 

In the following the atoms within the distorted 
zone in this sense are disregarded, and only the 
indirect effect of the distortion, through its 
elastic strain field, is taken into account. 

Consequently, we may assume that in Eq. (3) 
u, is small in comparison to the interatomic dis- 
tance. With the notation 


G=>., exp (27iR-r,), (4) 


the Taylor expansion yields: 


G= exp (27iR-a,) 
un exp (27iR-a,)) 


(5) 


The first term of (5) represents the ordinary 
scattering by a regular crystal and can be dis- 
regarded here. As we are interested only in the 
scattering outside the lines, the second term may 
be treated separately. 

In order to bring the second term of (5) into a 
more convenient form, we use a transformation 
formula given by Ewald.® If A, are the lattice 
vectors of the reciprocal lattice and v the volume 
of the unit cell, then 


1 
un exp U(R-A,), (6) 


where U(R) is the Fourier transform of a function 
*P. P. Ewald, Géttinger Nachrichten 55 (1938). 
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u(r) which interpolates the values u, at the lattice 
points: 


U(R)= f u(t) exp (7) 


and 


(8) 


This transformation may be considered as a 
three-dimensional generalization of Poisson’s sum 
formula. Patterson’s representation of the inter- 
ference function of perfect crystals® can be 
recognized as a special application of the gener- 
alized Poisson formula. 

The Eggs. (6)—(8) are evidently still applicable 
if the quantities u, are vectors, u,, if only the 
functions u(r) and U(R) are replaced by vector 
functions u(r) and U(R). Equations (6)-—(8) 
become : 


1 
> u, exp (27iR-a,) =— U(R—A,), (6a) 
n v : 


u(r =a,) =Uy. 


U(R)= f u(t) exp (Ta) 


u(r=a,) =Un. (8a) 


The Eqs. (7) and (8) do not determine uniquely 
the function U because any function which 
interpolates the values u, is acceptable for u(r). 

The disorder scattering function G; is built up 
of a number of identical functions having their 
respective origins at lattice points of the re- 
ciprocal lattice. The advantage of this transfor- 
mation is that the function U(R) often decreases 
so rapidly from R=0, that near a lattice point A; 
of the reciprocal lattice, all terms of thé sum (6) 
can be neglected except the term U(R—A),). 


As an example, the case of spherical dilatation 
will be considered. If one or several atoms are 
inserted into the lattice in addition to the regular 
atoms or if a larger atom is substituted for a 
regular one, then the neighboring atoms will be 
pushed radially outwards. In most metal crystals 
each atom is surrounded by a sizable number of 
equidistant atoms so that it is reasonable to 
replace the central perturbation by a spherical 
cavity with a pressure ~;. Furthermore, the 
resultant displacements will be calculated as if 


* A. L. Patterson, Phys. Rev. 56, 972 (1939). 
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the crystal were isotropic, because the results will 
be interesting mainly for their order of magni- 
tude. According to Love,!® the radial displace- 
ment u in a shell bounded by concentric spherical 
surfaces of radii ro and 7:(r9>71) is given by 


1 puri® 1 1 
(9) 
3A+2u 


4u ro* r? 


In our case the crystal size ro is large as compared 
to the radius of the distortion r; so that the first 
term in Eq. (9) can be neglected. Similarly, we 
disregard the term r;* in the denominator of the 
second term of (9). The displacement u; at the 
inner surface is 


ui =(ri/4u)pi, (10) 
and we get instead of (9) | 
u =u, (r1?/r?). (11) 


Inserting (11) into (7a), we obtain 
r 
U=u;r;? f exp (2riR-r)dv. (12) 
r 


Because of the spherical symmetry of the dis- 
placement, U must be parallel to R: 


cos (r, R 
f exp (27iR-r)dv, 
r 


(13) 


R \ 2mriR 2mroR 


This function vanishes at R=0, but beyond a 
maximum at R=1/2ro, the first term rapidly be- 
comes predominant. Now, the central spot of the 
interference funetion Gp which gives rise to the 
x-ray line itself, extends to about R=2/ro. As the 
disorder scattering cannot be measured in the 
immediate neighborhood of the line, we may limit 
ourselves to such points which are well outside 
the line. In the following, this domain will be 
called “near the lattice point” (of the reciprocal 
lattice) or, ‘“‘near the line.’’ With this restriction, 


Ur becomes 
sin 
Ur (14) 


0A. E. H. Love, Mathematical Theory of Elasticity 
(Cambridge University Press, New York, 1934), p. 142. 
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and, by Eqs. (6a) and (5) 


G,= cos (R, 
v 
sin 2rr,|R—A,| 
[R—A,|? 


Now, the distance between two neighboring 
lattice points of the reciprocal lattice is of the 
order of the inverse interatomic distance. As 7; is 
of the order of this latter quantity, it may be seen 
from (14) and (5) that the bulk of the function 
which constitutes the contribution of a particular 
lattice point A; will be contained in the space 
between A; and its immediate neighbors. Near 
any lattice point, we may disregard the contri- 
butions of all other lattice points, because of 
the rapid increase of the individual function 
U(R-—A,). 

The surrounding of the lattice point (000), i.e., 
of the direct transmitted beam, assumes a special 
place. By Eq. (5), the disorder scattering is pro- 
portional to R, and furthermore }-u, vanishes 
because the center of gravity is at rest, so that we 
should expect the disorder scattering to tend 
toward zero near the direct beam. The physical 
significance of this fact was explained by Ewald" 
by the observation that even an asphalt road 
makes a fairly good mirror when observed at a 
very oblique angle. Yet, according to Eq. (15), 
the predominant term of the sum is that which 
corresponds to A,=0; R becomes equal to 
(R—A,): 

2uir; sin 


’ 16 
, v R (16) 


and G;, tends toward a constant value at R=0. 
This is true, as explained above, only for points 
not too close to the line (000) itself, but it 
describes the experimental form of the disorder 
scattering curve. Indeed, Guinier’s experiments 
show that the background intensity of cold 
worked copper becomes constant near the scat- 
tering angle zero. As the Compton radiation as 
well as the temperature scattering tends toward 
zero at small angles, it is probable that the effect 
is caused by a perturbation similar to that con- 
sidered here. 


uP, P. Ewald, Handbuch der Physik (Berlin, 1927), 
Vol. 24, p. 271. 
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Equation (14) predicts a strong increase of 
intensity near a reciprocal lattice point. This, too, 
is in agreement with Guinier’s observations. The 
exact form of the intensity curve for a poly- 
crystalline aggregate would be given by the 
average of |G,|? over all crystal orientations, the 
absolute value |R| being fixed for a given 
scattering angle. But this integration would 
not have a result comparable to Guinier’s experi- 
ments. Indeed, G; is proportional to the cosine of 
the angle between R and (R—A,) in the example 
under consideration, and this is owing to the 
spherical symmetry of the perturbation assumed 
in this section. We cannot expect the local 
perturbation created by plastic deformation to 
have a spherical symmetry, so that our example 
would not be typical for the general case. 


3 


The results of the preceding section are sus- 
ceptible of considerable generalization. In this 
section, an infinitely extended body with a local 
perturbation of general character and no forces at 
infinity will be considered. The action of the local 
distortion will be replaced by a system of volume 
forces f(r), confined within a volume of atomic 
order of magnitude. The resultant force, 


f f(r)dv, 


must vanish if there are no boundary stresses. 
The laws of isotropic elasticity will again be used. 
This method differs from that of Section 2, 
where the distorted zone was replaced by a 
cavity, in that the atomic configuration within 
the distorted zone is replaced by one which is 
formally in agreement with elastic equilibrium 
and the fictitious volume force. However, because 
of the small number of these atoms, no appreci- 
able difference between the two methods should 
arise. 
If the displacement u(r) is represented by a 
Fourier integral 


f U(R) exp (—2iR-r)dV, (17) 


then the inversion of Eq. (17) is Eq. (6a). By 
differentiating under the integral sign, we obtain 


from Eq. (17): 
— te" f R°U exp (—2riR-r)dV, (18) 

grad div u 
f R-(R-U) exp (—2iR-r)dV. (19) 


The force function f(r) will also be represented by 
a Fourier integral : 


f(r) = f F(R) exp (2xiR-r)dV, (20) 
with the inversion: 


F(R) = f exp 


The equation of elastic equilibrium: 
wAu+(A+yz) grad div u+pf=0, 
becomes, by Eqs. (18)—(20) : 


(22) 


fi exp (—27iR-r) 
—42°R-(R-U)+pF]dV=0. (23) 


By virtue of a general theorem, the transform 
must vanish: 


(24) 


The solution of this vector equation is: 
p /RX(FXR) R-(F-R) 
+ ). @s) 
A+2p 


We are interested in the behavior of U at small 
values of |R|. Expansion of Eq. (21) yields: 


The first term is the resultant force, and there- 
fore, vanishes. Introducing the expansion (26) 
into Eq. (25), it can be seen that U behaves like 
1/R for small values of R. 

This is the generalization of the result of 
Section 2. In particular, it can be seen from 
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Eqs. (5) and (6a) that near the lattice point 
(000), the scattering amplitude and consequently 
the intensity tend toward a constant value, in 
agreement with Guinier’s observations. 

The results of this section are valid only for 
infinitely extended crystal grains. When the grain 
is finite, but large enough to give rise to sharp 
lines, the statements of this section remain valid 
only for points well outside the line, as was shown 
in Section 2. 


Near a lattice point A,, the scattering ampli- . 


tude will be given by 


Now, if the force system f has a privileged 
crystallographic direction, then we must expect 
to obtain different scattering near different 
lattice points of similar absolute distance | A,|, 
because of the direction cosine appearing in 
Eq. (27). Conversely, observations of the disorder 
scattering near different lines should be able 
to show the crystallographic symmetry of 
distortions. 

A disorder scattering of the _ under con- 
sideration must affect the measured line intensity 


and line width. However, conditions are more. 


complex within the line itself, because the grain 
size and the slowly variable stresses account for 
part of the line broadening.® 
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A study of the AlO spectrum indicates that a predissociation of the type one 6 takes place 
in the lower electronic state of the green system. It is suggested that the y- and e-bands of NO 
do not belong to a single system; the y system shows a strong perturbation of the v’ =4 level. 
An accidental predissociation of the vibrational type is found in the spectrum of Tez, and the 
analogy of the three molecules S2, Sez, and Tez is strengthened. 


AlO 


NLY a few examples of predissociation in- 
volving the lower state of a molecule are 
known, among which the best studied case is 
that of HgH. However, in the type III of pre- 
dissociation (by rotation) exhibited by HgH the 
widening of the predissociated levels is such that 
they cannot be observed. Holst! has reported 
predissociation effects in two systems of emission 
bands of AIH due to transitions ending on the 
unstable "II level. But, so far as we know, this 
example has not been studied in detail. The 
observation of the predissociated lines in the 
above-mentioned examples is difficult on high 
dispersion spectrograms, on account of the in- 
tensity drop in the centers of the lines. 
It is well known that, as a result of the rise in 
total absorption in the incompletely resolved 


1W. Holst, Thesis (Trondheim, 1936). 


predissociated lines (Loomis and Fuller),? pre- 
dissociation may, under certain conditions be 
observed in absorption even with low dispersion. 
This criterion has been applied previously to the 
investigation of induced predissociation. In the 
case of natural predissociation, the criterion was 
first applied to the Hz: spectrum by Beutler, 
Deubner, and Jiinger,? and to the Sz spectrum 
by Rosen and Neven.‘ The latter example has 
recently been investigated in detail by Herzberg 
and Mundie.° 

In emission, the widening of the lines of in- 
completely resolved bands manifests itself in 
different manners according to the pasition of 

2 F. Loomis and H. Fuller, Phys. Rev. 39, 180 (1932). 

*H. Beutler, A. Deubner, and H. O. Jiinger, Zeits. f. 
Physik 98, 181 (1935). 

*B. Rosen and L. Neven, Comptes rendus 203, 662 
(1936); J. Chim. Phys. 35, 58 (1938). 


G. Herzberg a % J Mundie, Phys. Rev. 8, 263 
(1940). 


SPECIAL CASES OF PREDISSOCIATION 


the first predissociated line relative to the head. 
If this first predissociated line is far removed 
from the head, a sudden breakdown in intensity 
should be observed, the band being thus shorter 
than the non-predissociated ones. If the pre- 
dissociation occurs at a J-value nearer to the 
head, the effect of intensity drop may be partly 
compensated by the superposition of widened 
neighboring lines. Finally, a predissociation oc- 
curring in the region of the head should result in 
an abnormal weakening of the whole band, 
which should then appear diffuse and degraded 
in both directions. As the J-value corresponding 
to the head varies from one band to another, the 
extension and intensity of the bands in a given v”’ 
progression varies in a characteristic manner 
according to v’. We have observed these charac- 
teristic changes in the spectrum of AlO and have 
applied them to the study of predissociation. In 
the v’’=7 progression, for instance, the following 
observations can be made. The (4,7) band is 
abnormally weak; it is diffuse and degraded in 
both directions. The (5, 7) band has a normal 
intensity in the region of the head, but this band 
is short and slightly degraded to the violet in 
addition to the normal degradation to the red. 
The (6, 7) and following bands have normal 
intensities and become gradually longer, the 
breakdown of intensity occurring at a greater 
distance from the head. Owing to superpositions, 
an accurate determination of the position of the 
breakdown is possible for the (9, 7) and (10, 7) 
bands only; it occurs in the former at a distance 
of 28 cm= from the head, in the latter at a 
distance of 37 cm~'. From the molecular con- 
stants of the green system of AIO it is found that 
the intensity drop corresponds to K = 44, v’’=7. 
In the same manner the observations in the 
vy’ =6 progression are explained by a predissocia- 
tion at K = 60, and those in the v’’ =8 progression 
by a predissociation at K=8, whereas all the 
bands of the v’=9 progression seem to be 
predissociated. 

The application of Herzberg’s theory to these 
results indicates that the predissociation takes 
place in the lower state of the 22 —* system of 
AlO, that it belongs to the type one b, and that 
the asymptote of the perturbing state is situated 
0.93 volt above the v’’=0, K=0 state of the 
lower *E level. It is probable that the predissocia- 
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tion observed by Coheur-Dehalu® in the (0, 0) 
band at K=128 belongs to the same limiting 
curve of predissociation. 

The lower state of the * system is generally 
considered as being the normal state of AIO, 
this system being observed in absorption in the 
laboratory under special conditions. On the other 
hand the vibrational levels of this state are 
known up to 2.2 ev.’ It follows that one of the 
atoms resulting from the dissociation of AlO in 
the normal state must be excited; probably the 
end products of the dissociation are Al(?P) and 
O('D). A detailed account will be published 
elsewhere. 


NO 


Herzberg and Mundie’ have suggested that the 
bands previously attributed to the y- and e- 
systems of NO belong to a single system, the e- 
bands corresponding to transitions to the levels 
v>4 of the *2 state. The considerably greater 
intensity of the e-bands is explained as a pre- 
dissociation effect. 

The observation by Guillery® of y-emission 
bands corresponding to v>4 was not considered 
by Herzberg and Mundie as a conclusive ob- 
jection to their suggestion. But Guillery’s obser- 
vation, together with several other facts, such as 
the general character of the absorption bands of 
NO studied by Leifson® made a re-examination 
of this question desirable. Measures by P. 
Migeotte and the writer of all the NO bands 
excited in emission have enabled us to improve 
the vibrational classification of the. y-bands 
previously published by Guillery on the basis of 
the measures by Bair and Deslandres. A con- 
siderable perturbation of the vibrational type 
(displacement of about 90 cm~") affects the v’ =4 
level, but the y-bands attributed by Guillery to 
v’>4 appear with great intensity. An elaborate 
discussion to be published soon proves conclu- 
sively that these bands belong to the y-system. 
In addition, our observations of the emission 
bands in the region of 1900A exclude the possi- 
bility of considering the e-bands as part of 


a a ) Coheur-Dehalu, Bull. Acad. Roy. Belgium 23, 604 
7E, Bodson and F. Coheur-Dehalu, Bull. Acad. Roy. 
Belgium 23, 408 (1937). 
M. Guillery, Zeits. f. Physik 42, 121 (1927). 
*S. W. Leifson, Astrophys. J. 63, 73 (1926). 


ileged 
xpect : 
ferent | 
| 
1g in 
order 
able 
y of 
| 
more. 
grain 
it for 

1945 | 

pre- 
s be | 
sion. 
> the 
1 the 
was 
itler, 
‘(rum 
has 
berg 
f in- 
f in 
n of 
2). 
its. 
, 263 

\ 


126 L. PAGE AND N. 


the y-system. Hence the hypothesis of Herzberg 
and Mundie and its conclusions should be 


abandoned.!°® 
Te. 


In collaboration with R. Migeotte we have 
recently studied the ultraviolet part of the Tes 
spectrum." The following observation concerning 
a predissociation effect is of interest. The ultra- 
violet end of the Tes spectrum consists of the 
v’’ =0 progression. Theoretically this progression 
should have a regular intensity distribution with 
a single maximum. This is what we observe in 
absorption, the maximum being found at v’ = 15. 
In emission, however, the bands v’=14 and 
v’=16 are abnormally weak (for the number- 
ing of the Tez bands).” It seems that we are 
confronted with a predissociation phenomenon 
which can be accounted for as an accidental 
predissociation of the vibrational type analogous 
to a vibrational perturbation affecting a whole 
band. 

10 We have recently published a detailed investigation of 
the various ultraviolet band systems of NO (Bull. Soc. 
Roy. Sci., Liége, Belgium, pp. 40 and 49, 1945). Subse- 
quently we have learned that three papers had been 
recently published on the NO molecule, one by Gerd, 
Schmid, and Szily (Physica 9, 144 (1944)); the others by 


Gaydon (Proc. Phys. Soc. 56, 95, and 160 (1944)). As far as 
predissociation is concerned, these authors reach the same 
conclusion as we do. 

_™R. Migeotte and B. Rosen, Bull. Soc. Roy. Sci., Liége 


(in press). 
-  #R, Migeotte, Bull. Soc. Roy. Sci., Liége 13, 48 (1942). 
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The unresolved isotopic bands v’=14 and 
v’=16 are not only weakened, but they present 
in emission a very pronounced asymmetry. In 
connection with the observations of Olsson," 
who states that the positions of the individual 
absorption bands are quite normal, this means 
that the weakening by accidental predissociation 
is different for different individual isotopic bands. 
No such asymmetry is observed in absorption. 

Two regions of predissociation are known in 
the Sz and See spectra. In Se, predissociation 
occurs at the levels v’>10 and v’>17° whereas 
in Sez it occurs at the levels v’>10 and v’>22.¥ 
In the spectrum of Tes only one predissociation 
region, setting in at v’=21, has been reported." 
Our new observation indicates that the Te, 
molecule also can predissociate at intermediate 
v’ values. However, contrary to the other mole- 
cules of the same group, the predissociation of 
Te: at intermediate v’ values is only possible 
through an intermediate state, by an accidental 
predissociation of the vibrational type. The close 
analogy existing between the spectra of the three 
molecules is thus strengthened, especially since 
considerable perturbations of the vibrational 
type are known for numerous excited levels of 
S2 and 


43 E. Olsson, Thesis (Stockholm, 1938) ; Zeits. f. Physik 95, 
215 (1935). 

4 B. Rosen, Physica 6, 205 (1939). 

18 EF. Hirschlaff, Zeits. f. Physik 75, 315 (1932). 
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New solutions of the space charge equation are obtained, which converge much more 
rapidly than Langmuir’s solution in the important case where the radius of the outer electrode 
is large compared with that of the inner electrode. 


HE theory of the limitation of a current 
by space charge in an evacuated region 
containing ions of one sign only was given first 
by Child! for parallel plane electrodes and later 


1C. D. Child, Phys. Rev. 32, 498 (1911). 


by Langmuir? and Langmuir and Blodgett? for 

the far more important case of coaxial cylindrical 

electrodes, as found in the usual construction of 
*T. La ir, Phys. Rev. 2, 450 (1913). 
ogmuir and K. B 


37. Langmuir . Blodgett, Phys. Rev. 22, 347 
(1923). 


SPACE CHARGE BETWEEN COAXIAL CYLINDERS 


the diode vacuum tube. In the first case the 
solution is obtained easily in closed form, but in 
the second a non-linear differential equation is 
encountered which can be solved only in series. 
Langmuir’s solution is a power series in the 
logarithm of the ratio of the distance r from the 
common axis of the electrodes to the radius a 
of the inner electrode or cathode. While this 
series converges satisfactorily for small values 
of r/a, it converges very poorly in the much 
_more important region where r/a is large, and 
special approximate methods had to be em- 
ployed for values of r/a greater than 66. The 
object of the present paper is to develop a 
solution which converges the better the larger 
is the ratio r/a. In fact, five terms of the solution 
obtained give an accuracy of two or three parts 
in ten thousand for all values of r/a greater 
than 3. 

If e/m is the ratio of charge to mass of the 
ions (electrons) in the tube, V(r) the excess of 
the potential over that of the cathode, p(r) the 
charge density, and 7; the current per unit legnth 
of the coaxial electrodes, the equations to be 
satisfied are, in Heaviside-Lorentz units, 


2e 
r= — 


(1) 


(2) 


(3) 
subject to the boundary conditions 
dV/dr=0, p=, (4) 


at r=a. The boundary conditions imply satura- 
tion current and negligible initial velocities. 
Eliminating p and 7, we find, if we put n=a/r, 


2e\'(—V)! 
rg(n) 
where g(n) satisfies the non-linear differential 
equation 


(S) 


(6) 
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The function g is identical‘ with Langmuir’s 
It is obvious that g=1 satisfies (6). Hence to 
obtain a solution which converges rapidly for 
large r or small n, it would be natural to try the 
power series 


However, substitution in the differential equa- 
tion shows that a power series solution in 7 
does not exist. Hence we change the independent 
variable from 7 to z, getting in place of (6) 


(= 


2 


Ws dg 
dz 
dz dey 


dn dz 
and assume a power series solution in z, to wit: 


Substituting in the differential equation and 
equating the coefficients of a; on both sides, we 
find that a, is arbitrary provided 


0, 
dn 
the solution of which is 
in?!+A 2??, 


2+i(2)! 2—1(2)! 


’ 


where 
p2= 


Now put 
2)! 


( 
log y+a. 


Then a power series solution in 
zi=cye*, 
exists, of the form 
g=14+- (21422) + (@1121? +02222") 
+ (@11121° +01 1221722 
(7) 


This is a complete solution, since it contains the 
two arbitrary constants c and a. 


‘The reason for this departure from an established 
notation is that our treatment leads to a differential 
penny 7 , and it would be awkward and unusual to 

dent variable as the square of a function. 
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To evaluate the coefficients in (7), it is con- 
venient to change the independent variable in 
(6) from 7 to y. Then 


“( »| Bal (8) 
dy 
and we find for the first five terms of the solution’ 


g=1+2cy cos 8 
2(2)! 


1 
sin 28+-| 
33 3 


5 
— cos 28+ 
33 


19(2)! 
sin 38 
2376 


7541(2)) 
—————_ si 
1,685,772. 
82(2)) 17 


sin vee, (9) 
16,929 11,286 


n 48— cos 28 
33,858 


When we attempt to determine the constants 
c and @ from the boundary conditions g=0 and 
_ dg/dy=0 at y=1, we find that this series does 
not converge rapidly enough for this value of y 
to ensure the desired accuracy of one part in a 
thousand from the five terms computed. Instead 
of calculating further terms we have developed 
a second series which converges very rapidly in 
the neighborhood of y=1, and then made use 
of the principle of analytic continuation to fit 
the series (9) to it. This second series is a power 
series in x=1—y=1—7?!, 

Changing the independent variable in (8) from 
y to x we have 


9 dg} 
O<x<1. (10) 
2 dx 


* After expressing the solution in trigonometrical form, 
the numerical coefficients were checked by substituting 
back in the differential equation. 
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In order to satisfy the boundary conditions at 
the cathode, the function g must be of the form 


9 


We find the following recursion formula for the 
coefficients: 


—2 ]by-1 
—[3p?—10.756,667p+10.41 ]b,_» 
+[(5b3—8b2+3b,)p 

— (6563 —68b2+16.5b1) 
+[(6b,—10b;+4b.)p 

— (102b,—125b3+36.5b2) 
+[(7bs;—12b4+5bs)p 

— 

— (200bs— 287bs+100.5b4) 


Joy (11) 


where first differences are constant in the coeffi- 
cient of p and second differences in the other 


TABLE I. 


g from (12) g from (14) 


0 

0.2793 
0.5171 
0.6671 
0.7666 
0.8363 
0.8870 
0.9254 
0.9550 


= 
» 


NK 


; 17 
——— cos 
594 
1 
+—— cos 8+—— sin 8 
396 264 
2386 
——— cos 48 
421,443 
f 
t 
p 
= x 
0 te 
0.37004 fy 
| 0.51925 0.48075 0.5170 
0.60315 0.39685 0.6670 Ci 
0.65800 0.34200 0.7667 
0.69715 0.30285 0.8363 g 
0.72672 0.27328 0.8870 
0.75000 0.25000 0.9254 
0.76888 0.23112 0.9549 
0.21544 0.9782 
0.20218 0.9970 
0.19079 1.0123 
0.13572 1.0716 
0.08550 1.0947 
0.05386 1.0846 
0.03393 1.0634 
0.02137 1.0422 
0.01347 1.0253 . 
1280 0.00848 1.0135 
2560 0.00534 1.0060 
5120 0.00337 1.0019 
10240 0.00212 0.9998 
cla 20480 0.00134 0.9991 
40960 0.00084 0.9990 
\ 


is at 
form 


the 


term. This gives for the solution 


9 
—0.200,000x —0.110,833x? 
4 


—0.065,303x* —0.042,006x* 
—0.028,876x' —0.020,858x*® 
—0.015,647x? —0.012,091x8 
—0.009,568x° —0.007,720x!° 
—0.006,330x!! —0.005,261x!? 
—0.00442x!8 —0.00375x™ 
—0.00321x%—---]. (12) 


This series is superior to that given by Langmuir 
in that the variable goes to unity instead of 
infinity as r/a increases without limit. 

After a few trials it appeared that x=2/3, 
y=1/3 was the best point at which to fit the 
complete solution (9) to the solution (12) satis- 
fying the boundary conditions. At this point 


dg 
g= 0.78220, — 0.59784, 
y 
which gives 
c= —0.9780, a=—22°.83, 


to one part in a thousand. Consequently, if we 


put 


6=93°.288 logio (1/y)+22°.83 (13) 


to avoid negative signs, the solution (9) satis- 
fying the specified boundary conditions at the 
cathode is 


g=1—[1.9560 cos 6]y+[0.63767 
+0.28985 cos 20—0.16396 sin 
—[0.00472 cos @—0.01002 sin 
— 0.05354 cos 39—0.02116 sin 36]y* 
+ [0.00276 — 0.00492 cos 26 
+0.01253 sin 20+0.01036 cos 48 


+0.01158 sin 40]y'—---. (14) 
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These five terms give an accuracy of better 
than one part in a thousand for all r/a greater 
than three. In fact, for r/a greater than ten the 
first four terms are sufficient, for r/a greater 
than fifty only three terms are needed, and for 
r/a greater than five hundred two terms suffice. 


1.00 


Fic. 1. Value of the function g. 


We have calculated g for a number of values 
of r/a from both the “‘near” formula (12) and 
the ‘‘far’’ formula (14) for the purpose of show- 
ing the excellence of the fit obtained. The results 
are included in Table I. It will be noticed that 
in the overlap extending from r/a=3 to r/a=9 
inclusive, where both series converge satisfac- 
torily, discrepancies are nowhere greater than 
two parts in ten thousand. Values of g calculated 
from (14) for larger values of r/a are also con- 
tained in the table. ; 

In Fig. 1 we have plotted the function g 
against the logarithm of 7/a, to illustrate how 
the function behaves in the range r=Sa to 
r=2560a. Beyond the latter point the curve 
oscillates about the straight line g=1 with ever 
decreasing amplitude. The function g becomes 
unity first at r=11.18a, again at r=9426a, and 
next at r=7,389,000a. In fact, g=1 for 


r=9426a[783.9]", (15) 


where n is any positive integer including zero. 
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A set of general equations of motion for classical dynamical systems of variable mass is 
derived. The variation of mass may be caused by either a continuous deformation and motion 
of the defining surface or a motion of the material points of the system, or both. Special cases of 


the general equations are also considered. 


PART 1 


N the discussion which follows, the term, 
dynamical system, S, will be used to denote 
the totality of those material points which lie 
within an arbitrarily chosen finite closed sur- 
face S. Neither the shape of the surface nor the 
volume V which it encloses need be independent 
of the time. The position vector q, of the centroid 
of the volume V relative to the origin 0 of an 
inertial reference triple (i, j, k) may also vary 
with the time. We shall, however, restrict our- 
selves to the discussion of systems such that, 
during a given time interval a<t<6, each point 
of the defining surface has a position vector 
which, together with its first two time deriva- 
tives, is a continuous function of the time. As a 
consequence of this restriction, the centroidal 
vector q, and its first two time derivatives will 
_also be continuous with respect to the time. 

We should note carefully the distinction be- 
tween the centroidal vector q, of the volume V 
and the centroidal vector q of the mass M which 
comprises the system. For the purpose of clarify- 
ing this distinction, there is represented in Fig. 1 
a system for which these two vectors are not 
equal. Here the surface ABCD comprises the 
surface S which defines the system. All of the 
mass M is contained within the shaded portion 
of the volume V. From the usual definitions of 
the centroids of masses and volumes, it follows, 
at once, that the centroidal vectors q, and q are 
unequal. The latter vector we shall refer to as the 
centroidal vector of the system. 

For the purposes of the present paper, it is 
convenient to divide dynamical systems into two 


* The opinions or assertions contained in this paper are 
the private ones of the writer and are not to be construed 
as official or reflecting the views of the Navy Department 


or the naval service at large. 


classes. The first of these comprises all those 
systems whose total mass remains constant; the 
second, those whose total mass varies with the 
time. Systems of the first class may be divided 
further into two subclasses. In the first subclass 
we place those systems whose mass, while con- 
stant in magnitude, does not retain its identity. 
In the second subclass, on the other hand, we 
place those systems whose mass not only remains 
constant but retains its identity as well. A 
system of each of these subclasses is illustrated 
in Figs. 2 and 3, respectively. In each of these 
figures, R represents a rigid uniform rectilinear 
distribution of matter whose velocity ve has a 
direction parallel to R itself. For both systems the 
defining surface S is rigid and, therefore, of 


Fic. 1. System in which the centroidal vector q, is not equal 
to the vector q. 
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SYSTEMS OF VARIABLE MASS 


constant shape. The surface S of Fig. 2 has a 
velocity vs which is equal to —vez. In accordance 
with our definition, the dynamical system repre- 
sented in this figure comprises all of the material 
points of R which are situated within the sur- 
face S. Obviously, as long as neither of the ends 
of R lies within S, the mass M of the system re- 
mains constant. The material points which com- 
prise the system, however, are changing con- 
tinually. This system belongs, therefore, to the 
first of our subclasses. The surface S which 
defines the system of Fig. 3 moves with a velocity 
vs which is equal to +vr. Consequently, the 
material points which comprise this system do 
not change with the time, i.e., the mass of the 
system retains its identity. We have, accordingly, 
a system of the second subclass. 

The translational equation of motion for sys- 
tems belonging to the second subclass is the well- 
known equation >>; F;= MA. In this equation the 
summation }>°; F; represents the sum of the ex- 
ternal forces, that is, the sum of those forces 
which arise from sources external to the defining 
surface S. The vector A is the acceleration vector 
for the centroid of the system. It is the purpose 
of the first portion of this paper to present the 
derivation of an analogous translational equation 
of motion of considerable generality for the class 
of dynamical systems whose mass is variable. 

We shall begin our derivation by considering 
three related systems S, S;, and S: whose de- 
scriptions follow. The system S is the system for 
which a translational equation of motion is to be 
derived. In the time interval a<t<b the mass MW 
of this system is assumed to be (a) finite and 
unequal to zero, (b) monotonic with respect to 
the time, and (c) a continuous function of the 
time together with its first and second time 
derivatives. Suppose now that the value of M 
corresponding to some arbitrarily chosen interior 
point ¢; of the interval a<t<b is Mz and suppose, 
further, that the material points comprising the 
system at this instant are identified so that their 
motion previous and subsequent to this instant 
may be studied. These material points comprise a 
second system S:; whose mass not only remains 
constant and equal to M; but retains its identity 
as well. For this system, then, the equation of 
motion is 


F2(t1; t) = Mo(ti) -Ao(ti; 2). (1) 
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In this equation and, indeed, throughout the 
ensuing discussion, the time ¢,; plays the role of a 
parameter. The symbol F,(¢;;¢) represents the 
‘sum of the external forces acting upon the system 
S: at the time ¢. To indicate that the various 


Fic. 2. System of constant mass, but the mass does not 
retain its identity. 


= 


Fic. 3. System of constant and identical mass. 


quantities which appear in this equation are 
those associated with the system S: we have used 
the subscript 2. For the two systems S and S; we 
shall make the assumption that at all points of 
the interval a<t<b except t=t, either (a), the 
material points of S: comprise a proper subset of 
the material points of S, or (6), the material 
points of S comprise a proper subset of the ma- 
terial points of S:. In accordance with this 
assumption, if ¢ is unequal to ¢;, one or the other 
of the relations, 


(M>M2), (2a) 
M(t) = M2(ti) — Milt; 2), (M<M), (2b) 


will be valid. In these relations the symbol 
Myi(t:;t) represents the mass of those material 
points which together with the material points 
of the system 5S: or the system S comprise, re- 
spectively, the material points of the system S 
or the system S:. These material points, of total 
mass M,, comprise the remaining one of the three 
systems S, S;,and S:. The situation corresponding 
to Eq. (2a) is represented graphically in Fig. (4a) ; 
that corresponding to (2b), in Fig. (4b). In each 
of these two figures, ABCDA, ABCEA, and 
AECDA represent the surfaces S, S:, and Ss, re- 
spectively. Taken in the same order these sur- 
faces define systems of mass M, M,, and M2. 
During the time interval a<t<b both M(?) 
and M,(t;) are continuous functions of the time / 
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Fic. 4a. A system corresponding to Eq. (2a). 


and have continuous first and second time de- 
rivatives. Accordingly, their difference M—M, 
has the same properties. We shall denote this 
function by the symbol ¢. By definition, then, 


(a<t<b). (3a) 


As a consequence of this definition and the 
Eqs. (2a) and (2b), we have the further relations 


¢'(t)=M'(t), (a<t<b); (3b) 
(a<t<b); (3c) 

o(ti; t) = Mi(ti; t), (M>M:2); (3d) 

(M=M:2); (3e) 
(M<M2). 

The primes in Eqs. (3b) and (3c) indicate differ- 
entiation with respect to the time 7. It will be 


noted that in these two equations the parameter 
t, does not appear. 

We consider next the centroidal vectors q(t), 
ri(t:; 2), and p(t:; ¢) of the systems S, and S2, 
respectively. These three vectors are defined by 
the relations 


i 
Fic. 4b. A system corresponding to Eq. (2b). 


_ The symbols dm and ¢ represent the differential 


of mass and its position vector, respectively, 
while the subscripts S, S,, and S2 indicate the 
systems for which the corresponding integrals are 
to be evaluated. As a consequence of the prop- 
erties of 1, M2 and the surfaces which define our 
dynamical systems, it follows that q and p are 
continuous functions of the time ¢ and possess 
continuous first and second time derivatives 
throughout the given time interval a<i<b. A 
similar statement is not true for the centroidal 
vector r:, for at the instant ¢; the mass M, is 
zero, and r, and its time derivatives are undefined. 
We shall restrict the properties of r, further by 
making the following assumptions: 


(a) The vector rors ri(t,;¢) exists. We shall de- 


note this saan by the symbol b(¢;). 


(b) The vector Lim r;’(t,; ¢) exists and is finite. 
tt) 


This vector we shall represent by the sym- 
bol w(t). 


(c) The vectors Lim r,’’(t,;¢) and 
t—+(t1 —0) 


Lim ry,” (t:;¢) exist and are finite. 
t—+(t1 +0) 


From the fact that the systems to which we 
restrict ourselves are bounded by finite closed 
surfaces, it follows that the vector b(¢,) is finite. 
At all points of the interval a<t<b except ¢; the 
vector r, and its first two time derivatives are 
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continuous functions of the time. This follows as 
a consequence of the properties of M, M2 and the 
surfaces which define our systems. 

Turning now to the various moments of mass, 
we find that in consequence of the properties of 
M, M2, q, and p, the two moments M2p and Mq 
are, together with their first two time derivatives, 
continuous functions of the time ¢ throughout the 
interval a<t<b. The moment M,r;, on the other 
hand, is discontinuous together with its first two 
time derivatives at the point ¢,;. These three 
moments are related according to the equations 


Mats) p(t; t) = -a(t) 

— 4), 
(M> M2); 
(M= M2); 


(Sa) 
- ptr; = (Sb) 
t) = M(t) -q() 
+Mi(ti; t)-ti(ti; 4), 

(M<M:2). (5c) 
Upon introducing the function ¢ which we defined 
above and a vector function r whose definition 
will follow, these three relations may be replaced 
by the single equivalent relation 


M2(ts) -p(tr; t) = M(t) -q(t) — 4), 
(a<t<b). 


The vector r is defined as follows: 
(a) For t¥t,, 


(6) 


r(t1; =ri(ti; 4); 


(b) For t=t;, 
r(ty ty) =b(t,) =Lim ri(ty t) 
tty 


r’ ty) =w/(t,) =Lim ri (t1; t) 

tt] 
r’’(ti; t1) is finite but otherwise arbitrary. From 
its definition it follows that the vector function r 
and its first time derivative are continuous at 
every point of the given time interval. As will be 
seen presently our derivation of the translational 
equation of motion for the system S does not 
require the function r’’ to be continuous at the 
point 

It is worth noting that although the function 
r’ may be discontinuous at the point ¢; the 
function (g-r)” is continuous at this point. This 
may be proved very easily. Thus, making use of 
the continuity of the functions ¢, ¢’, ¢’’, rand r’ 
at the point ¢;, we have at once 


tt) tt) 


tt} 
= (ty) - (ti) -w(ti) +Le(ti; -Lim (t1; 
tt) 


Since botli of the possible values of the vector Lim r,’’(¢,; ¢) are finite and since the scalar g(t; ¢,) is 


tt} 


equal to zero, it follows that the bracketed product on the right is a zero vector. Consequently we 


may write, simply, 


Lim -r(ti; = M"" (ty) - b(t) -w(h). 


(7) 


If now, we evaluate the function [¢-r]” at the point ¢; we obtain the equation 
= M"’ (ti) -b(t,) +2M'(t;) ty) (ty ty). 


Upon recalling, once more, that ¢(¢1; ¢:) equals zero and that r’’(¢;; ¢;) is finite, we see that the third 
term on the right side of this relation also vanishes. We may, therefore, write 


(8) 


The right sides of relations (7) and (8) are identical. Consequently, 


tt) 


(9) 


and the continuity of [¢-r]’’ at the point ¢; is established. 
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Differentiation of both sides of Eq. (6) twice with respect to the time yields the relation 
(10) 
which is valid throughout the interval a<t<b. The vector p” is, of course, the acceleration vector A, 
of the centroid of the system S:. In accordance with Eq. (1), therefore, we may replace the left side 
of (10) with the summation of the external forces acting upon the system Sb». With this substitution 
Eq. (10) becomes 
F,(t1; 4) (t)-q(t) +2M'(t) -q'(t) + M(t) 
Worthy of note. is the fact that whenever M> M, the summation F;,(t;; ¢) must include the forces of 
interaction between the systems S; and S:. In consequence of relation (8) the form taken by this 
equation when is 
1) =[M" (tr) (ts) (4s) + M(t) J 
(12) 
At this instant, moreover, the forces which act upon the system S: are identical with those which act 
upon the system S. In consequence we may replace the left side of (12) with the summation F(¢;) of 
the external forces acting upon the system S at the time ¢;. After making this substitution and deleting 
the subscript, 1, we obtain, finally, the equation 
F(t)=(M" (a<t<b). (13) 


This equation comprises one form of the translational differential equation of motion of the system S. 
Another form, in which the momentum L(t) of the system S appears explicitly, may be obtained very 
easily. If we differentiate Eq. (6) once with respect to the time we obtain the relation 


M2(ts) -p' (5) (45), (@<tso). (14) 


- The left side of this equation represents the linear momentum L2(¢,; ¢) of the system S, at the time ¢. 
For the momentum Lz, then, we have the formula 


Since =0 and r’(¢:; is a finite vector, it follows that the product t:) -r’(t:; represents 
a zero vector. For t=¢;, therefore, Eq. (15) takes the form 


Lo(ty; t1:) =[M’ (th) -q(t1) + M (tr) 
=(M' (tr) +M (tr) (4) (th) (a@<ti<d). (16) 


For this instant, moreover, the linear momenta L(t) and Le(t1; t1) are equal. In consequence we may 
replace the left side of (16) with L(¢,). This substitution, together with a deletion of the subscript, 1, 
yields for the momentum L(é) of the system S the formula 


L(t) (a<t<d). (17) 


The time rate of change of L obeys the relation 
(4) - bY) (a<t<b). (18) 
where the vector u(¢) is defined by the relation 
u(t) =b’(2). (19) 
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Upon comparing Eqs. (18) and (13) we see that the latter may be replaced by the equivalent relation 


F(t) (a@<t<d). 


(20) 


For many applications this form of the translational equation of motion will be found to be more 


convenient than the form (13). 


Strictly speaking, we have demonstrated the validity of Eqs. (13), (17), and (20) only for interior 
points of the interval a<#<b. If ¢, is an end point of the interval a more or less obvious modification 
of the above arguments leads to the same results. These equations are valid, therefore, throughout the 


interval a<t<b. 


PART 2 


We shall derive now an equation which con- 
nects the moment of the external forces relative 
to an arbitrary fixed point and the time deriva- 
tive of the angular momentum of the system 
relative to the same point. Let a fixed point J be 
chosen and let H(t), Hi(ti;¢), and 
denote, respectively, the angular momenta of the 
systems S, S,, and S; relative to this point. The 


latter systems are defined in the same manner as 


in the derivation of the translational differential 
equation of motion. At any instant in the interval 
a<t<b these momenta will be related by one of 
the following equations: 

H2(t1; 4) =H(t)—Hi(ti;t), (M>M:2); (21a) 

= H(t), (M=M:); (21b) 

4), (M<M2). (21c) 
We shall assume that H(¢) and H2(¢,; ¢) are con- 
tinuous functions of the time ¢, and that they 
possess continuous first time derivatives in the 
interval a<t<b. As a consequence of these as- 
sumptions the difference H—Hy, will have the 


same properties. Upon denoting this vector func- 
tion by the symbol h we may write immediately 


t)=H(t)—h(ti;t), (a@<t<b). (22a) 
From this defining relation and the relations 


(21a)-(21c), we may deduce readily the validity 
of the following equations: 


Hy’ (1; 4), (@<t<b); (22b) 
2), (M>M2); (22c) 
(M=M:2); (22d) 
h(t,; 4) = —Hi(t; 2), (M<M2). (22e) 


The primes in Eq. (22b) indicate differentiation 
with respect to the variable ¢. Since the system S: 
is one whose mass not only remains constant but 


retains its identity as well, the quantity Hy,’ 
satisfies the well-known relation 


t) = ¢). 


Here K, indicates the total moment, taken rela- 
tive to the fixed point J, of the external forces 
acting on the system S2. Upon eliminating H,’ be- 
tween this equation and Eq. (22b) we obtain the 
relation 


K,(t,; t) =H’ (t) —h'(t; 4), 


When t=?;, the forces acting upon the system S 
are identical with those acting upon the system 
S:. For this instant, therefore, the moments of 
these two sets of forces relative to the point J are 
equal, and Eq. (23) may be replaced with the 
relation 

K(t:) =H’ (ti) (a@<ti<b). (24) 
In this equation K(t,;) denotes the moment, 
relative to the point J, of the external forces 


acting on the system S at the time #;. The quan- 
tity 2(t:) is defined as follows: 


(ts) =h’ (th; =Lim 


(a<t<b). (23) 


if M> M, during the approach to the limit ; 


tt} 


if M < M; during the approach to the limit. 


If ¢; is an end point of the interval a<t<b, we 
can modify the above argument slightly and 
obtain again the results (24) and (25). Inasmuch 


as the parameter ¢; represents an arbitrary point 
of the interval a<t<b, we may drop the sub- 
scripts and write, more simply, 


K(t)=H'(t)—2(), (@St<b). (26) 


This equation comprises the second of the equa- 
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tions of motion which we wished to derive. We 
shall refer to it as the rotational equation of 
motion. 


PART 3 


The argument used in the derivation of Eq. (26) 
may be employed also in connection with the 
linear momenta L(¢), Li(¢:; and Le(t,; t) of the 
systems S, S; and S2, respectively. The result will 
be an alternative form of the translational equa- 
tion of motion. Corresponding to the Eqs. (21a)- 
(21c), we shall have now the equations: 


Lo(t1; =L(t)—Lilti; 4), (M>M2); (27a) 
Lo(t1; = L(t), (M=M:2); (27b) 
Lo(ti; 4), (M<M2). (27c) 


We shall assume that with respect to the time ¢ 
the momenta L and Lz are continuous functions 
and possess continuous first derivatives through- 
out the given interval. The difference of these 
two momenta comprises a function I(¢; #) which 
possesses the same properties. Corresponding to 
Eqs. (22a)—(22e), the equations now are: 


(@<t<b); (28a) 

4), (a<t<b); (28b) 

2) 4), (M>M2); (28c) 

(M=M2); (28d) 

(ti; =—Li(h; 4), (M<M2). (28e) 

Inasmuch as the momentum L» obeys the relation 
t) =Le'(ti; 4), 


where F; is the sum of the external forces acting 
upon the system S2, we may replace Eq. (28b) 
with the equation 


(@<t<b). (29) 


For the instant ¢; the sum of the forces acting 
upon the system S: will be equal to the sum of 
those acting upon the system S. Upon denoting 
the latter sum by the symbol F(t:) we may 
replace Eq. (29) for this instant with the equation 


F(t:) => L’(t) n(t;), (a <h <b). (30) 


The symbol n(¢) which appears in this relation is 
defined as follows: 


Fic. 5. One-dimensional distribution (R) of matter. 
n(t1) (t1; #1) =Lim [Ly'(t:; 4) 
tt) 


if M> M: during the approach to the limit ; 


n(t:) =I =Lim > (31) 
tt) 


if M < _M:2during the approach to the limit. ) 


The modified argument which must be employed 
when ¢; is an end point of the given interval 
(a<t<b) leads to the same results (30) and (31). 
Upon dropping the subscripts, as usual, we may 


write 
F(t)=L'()—n(t), (@<t<d). (32) 


This equation comprises an alternative form of 
the translational equation of motion (20). 


PART 4 


When used in connection with the kinetic 
energies T(t), T1(t:; t) and T2(t1; t) of the systems 
S, S;, and Ss, respectively, the same argument 
leads to the following equation: 


(@ShSb). (33) 


The symbol P(¢;) represents the power of all the 
forces, both internal and external, acting upon 
the system S at the instant ¢;. For the function 
u(ti), we have the following defining relations: 


=Lim (41; ) 


if M> M: during the approach to the limit. 


(34) 
=Lim 


if M <M; during the approach to the limit.) 
If we drop the subscripts from (33), we obtain the 


(31) 


(34) 


equality 


(@St<d). (35) 


In contrast with our previously derived Eqs. (13), 
(20), (26) and (32), this is a scalar rather than a 
vector equation. 


PART 5 


We shall derive next a useful one-dimensional 
form of the translational equation. In Fig. 5, R 
represents a one-dimensional distribution of 
matter whose linear density p is an analytic func- 
tion of the time ¢ and the distance ¢ from the ma- 
terial point E,,. The surfaces ABCDA, ABCEA, 
and AECDA comprise the surfaces S, S;, and So, 
respectively. We shall assume that during the 
interval a<t<6 all material points which pass 
through the surface S do so at the point B of this 
surface. In general, the velocity vz of the point B 
will not be equal to the velocity vs, of the 
material point B,, with which it coincides. The 
material point with which B coincides at the time 
t; is the point which we have referred to already 
as E,,. We choose the surface S2 so that its point 
E moves in coincidence with the particular ma- 
terial point £,,. This surface defines, therefore, a 
system whose mass remains constant and retains 
its identity, i.e., a system of the second subclass. 
The symbols Q, R and P denote, in this order, 
the centroids of the systems S, S;, and S:. For the 
remaining symbols in Fig. 5, the meanings are the 
same as in the previous portions of this paper. 
When t=?, the vectors OB(t), r(t,;¢) and (t:; ¢) 
are equal. Since, by definition, r(¢,; t:) =b(t:), we 
see that for this instant we may write 


b(t:) =OB(t,). 


Inasmuch as?, is an arbitrary point of the interval 
a<t<Q6, it follows that the vectors b and OB are 
equal to each other at every point of this interval. 
Accordingly, we may drop the subscript and 
write simply, 


b(t) =OB(Z). 
Differentiation of this equation with respect to 
the time gives for the vector u(t) =b’(¢) the result 
u(t) =b’(t) =va(t), (36) 
where Vz is the velocity of the point B. 


The vector w(#) may be obtained in the follow- 
ing manner. Let the position vector of the 
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centroid R relative to the point E be denoted by 
the symbol n(¢,; ¢). The vector r may then be ex- 
pressed in the form 


=OE(t; t)+n(t; 2). (37) 


If, now, the reference triple (i, j, k) has been 
chosen so that the vector j is parallel to R, the 
vector n may be replaced with the product m-j in 
which 1 is the j scalar component of n. With this 
substitution the equation for r becomes 


r(ti;t) t)+n(ti;t)-j. (38) 


Upon differentiating both sides of this relation 
with respect to the time ¢ we obtain the result 


r’(t,; 2) =OE’(t; 2) -j 
(39) 
where vz is the velocity of the point E and n’ -j is 
the velocity of the centroid R relative to E. In 
accordance with the definitions of w and r we 
may now write 
w(ty) (t1; 41) =Lim r’(t,; 2) 
tt) 
=Va(ti; n’ (ti; 
tt} 


Since, according to our present definition of the 
surface S:, the point E remains in coincidence 
with the material point £,,, it follows that the 
velocity vz of E is equal to the velocity vz,, of the 
point £,,. Furthermore, at the time ¢=4:, the 
material point £,, is identical with the material 
point B,,. For this instant, therefore, we have the 
result that vg=vs,,. With this substitution the 
above relation for w becomes 


W(t1) =VBm(t1) +{Lim n'(ti;t)]-j. (40) 


The scalar n obeys the defining relation 


0 
f dg, (41) 


in which 8 is the j scalar component of the vector 
EB from E to B. At any time ¢ the density p may 
be expanded in a power series of the form 
p(tr; & t)=p(ts; 0, t) +pr(ti; 0, 
po(t,; 0, t) 
all 


(42) 
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In this relation the subscripts 1, 2, 3, --- etc. placed upon the p’s indicate the first, second, 
third, --- etc. partial derivatives with respect to the variable ¢. Upon replacing p in relation (41) 
with its power series expansion and integrating the result term by term we arrive at the relation 
’ t) ’ 0, t) Bh; ’ 

2! 3 


+pi(t:; 0, t)- 3 4 


—n(t; t)- E (t1; 0, t)- B(ti; t) +p1(t1; 0, t)- 


ot; 0, 1): 
or, equally well, 


t)+pi(ts; 0, , 2) ] 


If we differentiate both sides of this equation with respect to the time ¢ we obtain the relation 
1; 4) t) 


t)+pi(ts; 0, t)- 0, d+ t): 


The limiting form of this maou as ¢ approaches ¢; and n and 8 approach zero is 
Lim Lim 2). (45) 


tty 
The limit on the right is simply the j scalar component - the velocity of the point B relative to the 
point E at the time ¢; when these points coincide. In consequence, we may write 
Lim n'(t,; t)=4 Lim B’ (ti; t) = —ve (tr; J, (46) 


tt] 


where vz and vB,, are the j scalar components of vz and vB, respectively. In view of this result, we may 
replace the Eq. (40) for w with the equivalent relation 


W(ti) =VBm(t1) (tr) J. (47) 
Dropping the subscript 1, as usual, we have for the vector w at any point of the time interval a</<d, 
w(t) (48) 


To obtain the desired one-dimensional form of the translational equation, we substitute into Eq. (20) 


the values of u and w given by relations (36) and (48) respectively. This substitution yields at once 


the relation 
F(t) (a<t<b). (49) 


Because of its importance we recall once more the fact that vz,, is the absolute velocity of the material 
points of R at the instant they cross the surface S. The forces F(t) must, of course, include the force of 
reaction at the point B. 


= . 
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PART 6 


When the system under consideration is one all material points of which have the same common 
velocity v(t), the translational equation of motion has an especially simple form. Let us assume that as 
jt, the mass M is greater than the mass M2. The quantity n which appears in (30) must then be 
calculated according to the first of the defining relations (31). We have 


n(t;) =Lim Ly’ 4) = Lim (ti; t)-v(t)+Mi(ts; 


tt) 


- Lim CMy'(t1; t)-v(t)]+0-v'(t). 
In consequence of the relation 


M(t) = M2(ti)+Mi(ts; t), (M> M2), 


we must have 
M'(t) => My (th t). 


In accordance with this result, Eq. (50) may now be written in the form 
= Lim 


The usual deletion of the subscripts gives us the more general relation 
n(t)=M’'(t)-v(t), (a<t<b). 
Upon substituting this result in Eq. (32) we have 
F(t) =L'(t) —n(t) (4) -v(t) M(t) -v'(t) ]— -v(0), 
F(t)=M(t)-v'(), (@<t<b). (54) 


This is the desired equation. Except for the fact that M is a function of the time, this relation is 
identical in form with the translational equation for systems of the second subclass. As usual, care 
must be taken to include in F(¢) the forces of reaction between the system and the mass which is 
entering or leaving the system. 


or, finally, 


PART 7 


In this, the final portion of this paper, we shall 
consider some of the factors involved in the 
choice of systems S, S,, and S: which satisfy the 
assumed conditions. 

Let us consider first the case presented by a 
system S whose change of mass throughout the 


time interval a</<b is caused by a flow of | 
matter through a single region o of its defining ' 


surface. In particular, we shall consider a system 
for which the flow of matter is outward, so that 
its mass M is a monotonic decreasing function of 
the time. Such a system is represented in 
Figs. (6a) and (6b). The surface ABCFDEA of 
Fig. (6a) represents the defining surface S for the 
time t=¢,. Since, when t=¢,, the material points 
which comprise the system S: are identically 
those which comprise the system 5S, we see that 
for this instant we may choose the defining sur- 


face S, to be identical with the surface S. In 
general, as the time ¢ increases, the surface S will 
change in both position and form. For the time 
t=t,+e, where ¢ is a positive infinitesimal satisfy- 
ing the inequality e<b—1,, its form and position 


A 


Fic. 6. System for which the flow of matter is outward 
through one region and the mass is a monotonic decreasing 
function of the time. 


cond, 
n 

(2a) 
(1) 
(45) 

= 
) the 
(46) 
may 

<i, 
(20) tt, 
FIG.6o FIG 6b 
49) 
rial 


are represented by the surface A1B,CiF,\D,EiA1 
of Fig. (6b). The region o of the surface S corre- 
sponding to the times ¢=¢, and t=¢,+¢ is repre- 
sented, respectively, by the surfaces ABC of 
Fig. (6a) and A,B,C, of Fig. (6b). Consider, now, 
the particular material points which, at the time 
t=t,, coincide with the portion ABC of the 
surface S. At the end of the time interval e these 
material points will, in general, occupy new 
positions. In their new locations in space they 
will define a surface which is represented in 
Fig. (6b) by the surface A’B’C’. This new surface 
need not, of course, be congruent with the surface 
ABC. Now an essential feature of any choice of 
the surface S: is that the material points which it 
includes at any time ¢ of the interval a<t<b 
must be identical with those which it includes at 
the particular time ¢=¢,. For the time t=t,+e, 
therefore, the surface S,; must be such that a 
portion of it is composed of the surface A’B’C’. 
With this restriction in mind we see that the 
surface A’B’C’F,\D,E,A’ of Fig. (6b) represents a 
suitable choice. In accordance with this choice 
the surface A’B’C’F,C,B,A,E,A’ may be taken 
to represent the surface S;. The remaining sym- 
bols in Figs. (6a) and (6b) have the same meanings 
as in the previous portions of this paper. 

The next case we shall consider is that of a 
system S whose change of mass throughout the 
interval a<t<b is caused by the simultaneous 
flow of matter through several regions of its 
defining surface. For the sake of simplicity we 
shall limit our discussion to a system for which 
(a), the number of such regions is two; (b), the 
flow of matter through both regions is outward; 
and (c), the regions are non-contiguous. The sur- 
face ABCFKJGDHINEA of Fig. (7a) represents 
the defining surface S of such a system for the 
time t=¢,. Mass is leaving the system through the 
non-contiguous regions ABC and GDH. At a 
later time, t:+e, the system S will have changed, 
in general, both its position and its form. 
The surface A,B, C\ Fi Ki, of 
Fig. (7b) represents the surface S in its new posi- 
tion and form. In this figure the regions through 
which mass is leaving the system are labelled 
A,B,C, and G,D,H,. To obtain the translational 
equation of motion for a system of this type we 
may proceed as follows: Let a surface S, be 
chosen subject to the condition that throughout 
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Fic. 7. System for which the flow of matter is outward 
through two regions and the mass is a monotonic decreasing 
function of the time. 


the interval a<t<b it must divide the system S$ 


into two systems S, and S each of which is of 
the type previously considered. This restriction 
will be met if we require the set of those material 
points which lie in the surface S, at any time of 
the interval a<t<6 to be identical with the set 
of those particular points which lie in this surface 
at the instant ¢;. In Fig. (7a) the surface KZUYN 
represents such a surface S,. At the end of a time 
interval e, the particular material points which 
coincide with the region ZU Y of this surface at 
the instant ¢; will have moved to new positions. 
The surface which they define in their new posi- 
tions is represented by the surface Z’U’Y’ of 
Fig. (7b). If our restriction upon S, is to be met, 
it is necessary that at the time t=,;+.e a portion 
of S, should be composed of the surface Z’U’Y’. 
Accordingly, a suitable choice of S, for this 
instant is represented by the surface K,Z’U’ Y’M,. 
In order to distinguish between the various 
scalars and vectors which are associated with the 
systems S, and S, we shall use the subscripts a 
and b, respectively. Thus, for example, the sys- 
tem of constant mass which becomes identical 
with the system S, at the instant t=¢, will be 
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SYSTEMS OF VARIABLE MASS 


denoted by the symbol S24, while the corre- 
sponding system for the system S, will be denoted 
by the symbol For the time 
and So. are represented by the single surface 
ABCFKZUYNEA; S, and S2,, by the single sur- 
face KJGDHINYUZK. The representation of 
the surfaces Sa, Si,a, S2,a, So, and for the 
time ‘=t:+¢€ appears in Fig. (7b): S, is repre- 
sented by the surface A,B,C, F,K,Z'U' Y'N,E\A 


Sre A'B'C'F,C,B 
A'B'C'F,K,Z'U'Y'N EA’, 
Ki 
Sur 


and 
Ki 
Corresponding to the systems S, and S, there 
will be two equations, each of the form (20). 
Thus, we shall have 


F,(¢) = La’ (¢) +[ua(t) — 2wa(t) ]- 
(a<t<b); 
F,(t) = Le’ (t) + [us (t) — 2wa(t) 


(55a) 


(a<t<b). (55b) 


We note now that in the summation F,(¢) of the 
external forces acting upon the system S, there 
will appear the forces of reaction across the sur- 
face S,. Reaction forces of like magnitude and 
direction but with opposite sense will appear in 
the summation F,(#). Consequently, if we add 
together the corresponding sides of Eqs. (55a) 
and (55b), we shall find that the reaction forces 
do not appear in the resulting equation. For the 
various scalars and vectors we shall have the 
relations 


F,(t) +F,(t) = F(¢), the sum of the external 
forces on the system S; 


=L’(t), the time derivative 
of the linear momentum of S; 


(56a) 


(S6b) 
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M,'(t)+M,'(t)=M'(t), the time deriva- 
tive of the mass of S; 


u,(t) -M,'(t) (2) - My’ (t) =u(t)- M(t); 
and 
Wa(t) - Ma'(#) +wa(t) (t)=w(t)-M’(t). (S6e) 


The relations (56d) and (56e) are to be regarded 
as defining relations for the vectors u(#) and w(?) 
in this case. In view of these relations, the 
addition of the Eqs. (55a) and (55b) will lead to 
the result 


F(t) —2w(t)]-M’(d), 
(a<t<b). (57) 


This, the translational equation of motion for 
this case, is in precisely the same form as the 
relation (20). 

An analogous argument will lead to a rotational 
equation of motion of the form (26). Thus, we 
shall have 


(a@<t<b), 
where, in this case, 


K(t) =K.(4)+K,(é), the moment of the 
external forces on S about the fixed 
point J; 

H’ (t) =H.’ +H,’ the time derivative 
of the angular momentum of the sys- 
tem S relative to the point J; 


(56c) 
(56d) 


(58) 


(59a) 


(S9b) 


and 
=da(t) (S9c) 


The relation (59c) is to be considered as a defining 
relation for the vector 4 in this case. Systems for 
which the flow of matter is inward through some 
regions of the defining surface S and, simul- 
taneously, outward through other regions may be 
treated in a similar manner. 

I am indebted to Dr. H. T. Muhly for several 
helpful discussions; to Lt. E. E. Betz, and Dr. A. 
E. Currier for a critical reading of my manuscript. 


|_| 
mS 
s of = 
tion 
e of 
set 
face 
‘ime 
hich 
e at 
ons. 
’ of 
tion 
'N,. 
ious 
the 
ts a 
sys- 
ical 
| be 


PHYSICAL REVIEW 


Letters to the Editor 


ROMPT publication of brief reports of important dis- 
coveries in physics may be secured by addressing them 
to this department. The closing date for this department is the 
third of the month. Because of the late closing date for the sec- 
tion no proof can be shown to authors. The Board of Editors 
does not hold itself responsible for the opinions expressed by 
the correspondents. Communications should not in general ex- 
ceed 600 words in length. 


The Theory of Diffraction . 


F. B. Prppuck 
Corpus Christi College, Oxford, England 
July 28, 1945 


HE use of Kirchhoff’s surface integral in the theory of 
diffraction has often been felt to be unsatisfactory, 
and many attempts at improvement have been made. The 
present state of the theory is described with references to 
Stratton.! These attempts, at the best, repair the defi- 
ciencies of the integral when the screen is perfectly absorb- 
ing or perfectly reflecting, and do not, I think, touch on 
the problem of a screen of any composition. Consideration 
of currents in sheets has led me to a more drastic process, 
valid when the wave-length is small. 

Suppose that a wave falls obliquely on a sheet of small 
thickness d, and of dielectric constant ¢ and conductivity y. 
If the wave-length is small, we can first solve the problem 
for an infinite sheet of thickness d and a plane wave, and 
then take the electric moments in the sheet as if they were 
finite and let them radiate like Hertzian oscillators. If 
¥ is very great, the moments are all in a thin layer on the 
surface. In that case, a primary wave A cos (wt—kro)/ro 
from Po (of any polarization) gives rise to a secondary 
wave A cos ¢ sin (wt—kr—kro)dS/drio at Q, where dS is 
an element of the diffracting screen at P, ¢ the angle of 
incidence, Io= PoP, and r= PQ. This is Kirchhoff’s formula 
in the most important case of direct vision, the only case 
in which it is confirmed. Consideration of space compels 
me to forego details, which will be found in a book on 
Currents in Aerials and High Frequency Networks to be 
published shortly by the Clarendon Press, Oxford. If we 


have a small hole instead of a screen, and a plane wave,’ 


the current per meter in the perforated sheet is the same 
as in the hole. The electrical effect of the two is nearly that 
of a complete sheet, which below the sheet is equal and 
opposite to that of the incident light. This is Babinet’s 
principle, a natural consequence of the smallness of the 
wave-length. I have investigated the diffraction of a sheet 
of dielectric constant « and conductivity y at any angle, 
but the results are complicated. In treating the problem 
this way we are treating the physical phenomenon as it is, 
and determining the electrical state of a diffracting screen 
at all points of it. Such a spoint of view must, I think, 
’ ultimately replace the delicate and restricted methods 
based on Kirchhoff’s integral. 


1Stratton, Electromagnetic Theory (McGraw-Hill Book Company, 
Inc., New York, 1941) pp. 460-70. 
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Ho-Distribution of Fission Fragments 


N. O. LAssEN 
Institute of Theoretical Physics, University of Copenhagen, 
Copenhagen, Denmark 
June 16, 1945 


S was reported earlier, the total charge of fission 
fragments has been determined by deflecting the 
fragments from a thin layer of uranium in the magnetic 
field of the Copenhagen cyclotron. The uranium was placed 
inside the cyclotron about 6 mm behind an internal target 
of beryllium which was bombarded by 6 Mev deuterons, 
thus being a strong source of neutrons. The dimensions 
of the uranium layer were: width 1.0 mm, length 20 mm, 
thickness 0.35 mg/cm*; it was placed parallel to the lines 
of the magnetic force and the fragments used were emitted 
nearly normally from the surface. They passed through a 
slit S; of variable width placed 10 cm from the uranium 
layer and through a second slit S; 10 cm further away; the 
latter had a fixed width of 3.0 mm and was covered by a 
foil of thickness 0.79 mg/cm? of mica and 0.25 mg/cm? of 
aluminium through which the particles entered into an 
ionization chamber filled with pure argon. The collecting 
electrode of the chamber was connected to a linear am- 
plifier and this in turn to a cathode-ray oscillograph, 
which was photographed on a moving film. The slit S, 
could be moved in a direction normal to the plane deter- 
mined by the uranium layer and the window S:2; the 
various positions of this slit then allowed particles with 
various values of Hp to pass from the uranium layer to the 
ionization chamber. For all positions of the slit was seen a 
background consisting of a great number of pulses due to 
recoiling argon atoms in the ionization chamber; but 
only when the slit S; was displaced between 7 and 14 mm 
from the plane through the layer and the mica window S, 
pulses of much higher energy corresponding to fission 
fragments occurred. 
As more fully reported,' the Hp-distribution curve for 


number of fragments 


Fic. 1. Energy distribution of fission fragments having — 
traversed 0.79 mg/cm? of mica. 


Ea se 
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the fission fragments was determined in this way. Later 
on, small improvements of the apparatus have made it 
possible to measure the Hp-distribution for each of the two 
groups of fission fragments separately. 

Figure 1 shows the energy distribution as obtained when 


Hox 10°5 oerstedxem 


Fic. 2. Hp-distribution of fission fragments. Full circles: light group. 
Open circles: heavy group. Crosses: Both groups. _ 


the slit S: was made so wide (7 mm) that the fission frag- 
ments entering through S, could have all possible values of 
Hp. The figure clearly shows that, apart from a very few 
exceptions, the pulses of energies higher than 34 Mev are 
caused by fragments of the light, most energetic group, 
while the fragments of the heavy group correspond to 
pulses between 9 and 34 Mev; pulses of energies lower than 
10 Mev are caused by coincidences between recoiling argon 
atoms. In the Hp-measurements the slit S,; had a width of 
1.0 mm and it was moved in steps of 1.0 mm; for the 
different positions of S, the number of fragments were 
counted, of course for equal neutron doses, and by means 
of the photographic record divided into two groups corre- 
sponding to the sizes of the pulses. In this way we get the 
number of fragments of each group as a function of the 
displacement of S; and we then calculate by means of 
the known magnetic field (which is not homogeneous) 
the Hp-distribution curves shown on Fig. 2. 

As seen, the two groups have nearly the same values of 
Hp, but the light group has evidently slightly higher Hp- 
values than the heavy one. Putting the most frequent 
value of mE to 8.5X10* mass units X Mev! we get for the 
total charge of the light and the heavy group 20¢ and 22., 
respectively, when ¢ denotes the electronic charge. This is 
in agreement with calculations by Professor Bohr* which 
predict an effective charge of fission fragments of about 20e; 
yet the agreement is not perfect, as the theory involves a 
higher effective charge of the light group than of the heavy 
one in contradiction to the present experiments. Neverthe- 
less, it is no real disagreement; as pointed out by Professor 
Bohr in the paper, the total fragment charge need not be 
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identical with the charge effective in electronic interactions. 

The author wishes to express his heartiest thanks to the 
Director of the Institute of Theoretical Physics, Professor 
Niels Bohr, and to Professor J. C. Jacobsen for their 
interest in the work and their great encouragement. 

1N. O. Lassen, Kgl. Danske Vid. Sels. Math.-Fys. Medd. 22, No. 2 
press. 

? Niels Bohr, Phys. Rev. 59, 270 (1941). 


The Synchrotron—A Proposed High 
Energy Particle Accelerator 
Epwin M. McMILLAN 
University of California, Berkeley, California 
September, 5, 1945 

NE of the most successful methods for accelerating 
charged particles to very high energies involves the 
repeated application of an oscillating electric field, as in 
the cyclotron. If a very large number of individual accelera- 
tions is required, there may be difficulty in keeping the 
particles in step with the electric field. In the case of the 
cyclotron this difficulty appears when the relativistic mass 
change causes an appreciable variation in the angular 

velocity of the particles. 

The device proposed here makes use of a ‘“‘phase stabil- 
ity’’ possessed by certain orbits in a cyclotron. Consider, 
for example, a particle whose energy is such that its angular 
velocity is just right to match the frequency of the electric 
field. This will be called the equilibrium energy. Suppose 
further that the particle crosses the accelerating gaps just 
as the electric field passes through zero, changing in such a 
sense that an earlier arrival of the particle would result in 
an acceleration. This orbit is obviously stationary. To show 
that it is stable, suppose that a displacement in phase is 
made such that the particle arrives at the gaps too early. 
It is then accelerated; the increase in energy causes a 
decrease in angular velocity, which makes the time of 
arrival tend to become later. A similar argument shows 
that a change of energy from the equilibrium value tends 
to correct itself. These displaced orbits will continue to 
oscillate, with both phase and energy varying about their 
equilibrium values. 

In order to accelerate the particles it is now necessary 
to change the value of the equilibrium energy, which can 
be done by varying either the magnetic field or the fre- 
quency. While the equilibrium energy is changing, the 
phase of the motion will shift ahead just enough to provide 
the necessary accelerating force; the similarity of this 
behavior to that of a synchronous motor suggested the 
name of the device. 

The equations describing the phase and energy varia- 
tions have been derived by taking into account time varia- 
tion of both magnetic field and frequency, acceleration by 
the “betatron effect” (rate of change of flux), variation of 
the latter with orbit radius during the oscillations, and 
energy losses by ionization or radiation. It was assumed 
that the period of the phase oscillations is long compared 
to the period of orbital motion. The charge was taken to 
be one electronic charge. Equation (1) defines the equilib- 
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rium energy; (2) gives the instantaneous energy in terms 
of the equilibrium value and the phase variation, and (3) is © 


the “equation of motion” for the phase. Equation (4) 
determines the radius of the orbit. 


Eo= (300cH)/(2zf), (1) 
E=£E,[1— (do) /(d@)], (2) 
d/,d¢ 
+Vsing 
300 dFy Eo df 
+ 
R= 300H. (4) 
The symbols are: 


E, =rest energy, 

V =energy gain per turn from electric field, at most favorable phase for 
acceleration, 

L =loss of energy per turn from ionization and radiation, 

H =magnetic field at orbit, 

Fo =magnetic flux equilibrium orbit, 

@ =phase of particle (angular position with respect to gap when elec- 


0), 
light velocity, 

R =radius of orbit. 

(Energies are in electron volts, magnetic quantities in 
e.m.u., angles in radians, other quantities in c.g.s. units.) 

Equation (3) is seen to be identical with the equation of 
motion of a pendulum of unrestricted amplitude, the terms 
on the right representing a constant torque and a damping 
force. The phase variation is, therefore, oscillatory so long 
as the amplitude is not too great, the allowable amplitude 
being +2 when the first bracket on the right is zero, and 
vanishing when that bracket is equal to V. According to 
the adiabatic theorem, the amplitude will diminish as the 
inverse fourth root of Eo, since Eo occupies the role of a 
slowly varying mass in the first term of the equation; if the 
_ frequency is diminished, the last term on the right fur- 
nishes additional damping. 

The application of the method will depend on the type 
of particles to be accelerated, since the initial energy will 
in any case be near the rest energy. In the case of electrons, 
E> will vary during the acceleration by a large factor. It 
is not practical at present to vary the frequency by such a 
large factor, so one would choose to vary H, which has the 
additional advantage that the orbit approaches a constant 
radius. In the case of heavy particles Eo will vary much 
less; for example, in the acceleration of protons to 300 
Mev it changes by 30 percent. Thus it may be practical 
to vary the frequency for heavy particle acceleration. 

A possible design for a 300 Mev electron accelerator is 
outlined below: 


H =10,000 
radius of orbit "= 100 cm, 
megacycles/sec., 


injection =300 kv, 
initial orbit =78 cm. 


Since the radius expands 22 cm during the acceleration, 
the magnetic field needs to cover only a ring of this width, 
with of course some additional width to shape the field 
properly. The field should decrease with radius slightly in 
order to give radial and axial stability to the orbits. The 
total magnetic flux is about ¢ of what would be needed to 


satisfy the betatron flux condition for the same final energy, 

The voltage needed on the accelerating electrodes 
depends on the rate of change of the magnetic field. 
If the magnet is excited at 60 cycles, the peak value of 
(1/f)(dEo/dt) is 2300 volts. (The betatron term containing 
dF,/dt is about } of this and will be neglected.) If we let 
V=10,000 volts, the greatest phase shift will be 13°. The 
number of turns per phase oscillation will vary from 22 to 
440 during the acceleration. The relative variation of E, 
during one period of the phase oscillation will be 6.3 per- 
cent at the time of injection, and will then diminish, 
Therefore, the assumptions of slow variation during a 
period used in deriving the equations are valid. The energy 
loss by radiation is discussed in the letter following this, 
and is shown not to be serious in the above case. 

The application to heavy particles will not be discussed 
in detail, but it seems probable that the best method will 
be the variation of frequency. Since this variation does not 
have to be extremely rapid, it could be accomplished by 
means of motor-driven mechanical turning devices. 

The synchrotron offers the possibility of reaching energies 
in the billion-volt range with either electrons or heavy 
particles; in the former case, it will accomplish this end at 
a smaller cost in materials and power than the betatron; 
in the latter, it lacks the relativistic energy limit of the 
cyclotron. 

Construction of a 300-Mev electron accelerator using the 
above principle at the Radiation Laboratory of the Uni- 
versity of California at Berkeley is now being planned. 


Radiation from a Group of Electrons 
Moving in a Circular Orbit 
Epwin M. McMILLAN 


University of California, Berkeley, California 
September 9, 1945 


A SINGLE electron of total energy E (rest energy = E,) 
moving in a circle of radius R, radiates energy at the 


rate L (electron volts per turn), given by: 
L=400r(e/R)(E/E,)', (1) 


where ¢ is the electronic charge in e.s.u., and E>>E,. In 
the synchrotron one has the case of a rather concentrated 
group of electrons moving in the orbit, and the total 
amount of radiation depends on the coherence between the 
waves emitted by the individual electrons. For example, if 
there were complete coherence, the radiation per electron 
would be N times that given by (1), where N is the number 
of electrons in the group. 

It is apparent from the above that an answer to the 
coherence problem is very important for any device in 
which groups of electrons are made to move in a circle 
with high velocity. This answer is given by a formula due 
to J. Schwinger (communicated to the author by I. I. 
Rabi). Schwinger’s formula gives the radiation in each 
harmonic of the period of revolution, in a form that allows 
easy computation for any distribution of electrons around 
the orbit. It leads to the following conclusions: 
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(a) Most of the energy in (1) lies in very high har- 
monics. 

(b) The coherence between the high harmonics from 
different electrons tends to become very small if 
the group has an appreciable angular spread. 

(c) The low harmonics are partially coherent, and 
give an energy loss per electron per turn (L’) 
depending on N, but not on E if E> >£,. 

(d) Because of fluctuations from a uniform distribu- 
tion, each electron also radiates the same amount 
L that it would if alone in the orbit. The total 
radiation per electron is thus L+L’. 


Values of L’ have been computed numerically from 
Schwinger’s formula for the case of N electrons covering 
uniformly an arc with an angular extent which is 1/m of a 
circle. This was done for m=2, 4, and 6; also the asymp- 
totic form for large m was obtained. These values can all 
be fitted within a few percent by the formula: 


L'~400x(e/R)X 2.4(m*8—1)N. (2) 


Applying (1) and (2) to the case where R=100 cm, 
E/E,=600, N=10" (1/60 microcoulomb, giving 1 micro- 
ampere at a 60 cycle repetition rate), and m=6, we get: 


L=780 volts, L’ = 1400 volts. 


Thus the radiation loss will not seriously affect the opera- 


tion of the synchrotron. Furthermore, L. I. Schiff has 
shown that the coherent part L’, which is mostly in the 
very low harmonics, can be strongly reduced by shielding. 


A Note on the Separation of Gases by Diffusion 
into a Fast-Streaming Vapor 
F. A. SCHWERTZ 
Mellon Institute of Industrial Research, Pittsburgh, Pennsylvania* 
December 8, 1944 
HILE it is well known that gas mixtures may be 
separated by diffusion into a fast-streaming vapor, 
the potential speed of the process is not generally appreci- 
ated. The simple calculation given below may serve to 
correct this situation. The calculation is meant to be 
illustrative and not quantitative in the sense that it would 
agree with experimental data. 

In Fig. 1, A and B are two channels connected by a cross- 
duct, C, of very small cross-sectional area. A pure gas is 
assumed to be traveling in A at a high velocity in the 
downward direction. Similarly, a vapor is flowing in B. The 
pressure differential across C is assumed to be such that the 
fluid velocity in the x direction is constant. To express the 
rate of transfer of the gas the equation: 


i=pv—D grad p (1) 


may be used.' In this equation i is the rate of gas transfer 
per unit area; p, the density of the gas; v, the fluid velocity; 
and D, the mutual coefficient of diffusion of gas and vapor. 
For the steady state div i=0, so that 


div (pv)—D div grad p=0. (2) 
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Taste I. Rate of diffusion, i, of Hz into water vapor. (i is expressed in 
cm?/cm? at 150°C and 1 atmos.) 


v (cm/sec.) 6 =0.01 cm 6 =0.1 cm 6 =1.0 cm 
—1000 2.34 000 
—100 120 0.234 000 
—10 156 12.0 0.023 
-1 164 15.6 1.20 

165 16.5 1.65 
+1 166 17.6 2.20 
+10 177 22.0 10.0 
+1 220 100 100 
+1000 1002 1000 1000 


TABLE II. Separation currents and compositions. (i is expressed in 
cm?/cm? Xsec. at 150°C and 1 atmos.) 


=0.01 cm b=0.1 cm =1.0 cm 

cm/sec. i He i i fue 
—1000 1.17 1.00 000 _ 000 _ 
—100 62 0.964 0.117 1.00 000 
—10 91.5 0.852 6.23 0.964 0.012 1.00 
-1 97.8 0.839 9.15 0.852 0.623 0.963 

0 98.5 0.837 9.85 0.837 0.985 0.837 
+1 99.3 0.836 107 0.826 1.63 0.678 
+10 107 0.827 16.3 0.678 10.0 0.500 
+100 162 0.678 100 0.500 100 0.500 
+1000 1001 0.500 1000 0.500 1000 0.500 


If p is considered a function of x only, and v,=v is con- 


sidered constant, then 


(3) 


If it is further assumed that the velocities in A and B are 
so large that the concentration of gas in B and of the vapor 
in A are substantially zero, (3) must be solved subject to 
the boundary conditions: 


p(x)=po at x=0, 


and 
p(x)=0 at x=b. 
The result is 
(4) 
For the current (1) gives 
(S) 
The limit of i as v approaches zero is ; 
io=poD/b. (6) 


Fic. 1. Gas and vapor streaming by a connecting duct. 
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In Table I the calculated values of i for the diffusion of 
hydrogen into water vapor are presented for D=1.65 
cm?/sec. at 150°C. 

The value of i for b=0.01 cm and v=0 corresponds to 
19,500 cu. ft. of hydrogen per square foot per hour. 

If now it is assumed that two gases will diffuse inde- 
pendently of one another, the ratio of the currents will be 


_ po2L1—exp (—vb/D1)] 


i po,iLl1—exp (—vb/Dz2) ] (7) 


and the fraction, f;, of the first gas in the separated portion 
will be 


fin {14 (= vb/D1)]\~ 


po,i1L1—exp (—vb/D2) 


THE EDITOR 


In Table II the values of i and f, =fHe are given for a gas 
mixture originally containing 50 percent Hz and 50 percent 
CO. D=0.32 cm*/sec. was used for the mutual coefficient 
of diffusion of CO, and water vapor. 

Table II shows that, at a fixed value of v, an increase in 
H: concentration may be effected merely by increasing the 
length, b, of the cross-duct, but only at the cost of a severe 
reduction in the gas current. 

The calculations suggest that a thin, perforated metal 
sheet would serve as the nucleus of an effective gas 


separator.” 
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